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Abstract 

We perform a comprehensive study of charged lepton flavour violation in Randall- 
Sundrum (RS) models in a fully 5D quantum-field-theoretical framework. We 
consider the RS model with minimal field content and a “custodially protected” 
extension as well as three implementations of the IR-brane localized Higgs field, 
including the non-decoupling effect of the KK excitations of a narrow bulk Higgs. 
Our calculation provides the first complete result for the flavour-violating electro¬ 
magnetic dipole operator in Randall-Sundrum models. It contains three contri¬ 
butions with different dependence on the magnitude of the anarchic 5D Yukawa 
matrix, which can all be important in certain parameter regions. We study the 
typical range for the branching fractions of p —> e'y, p —y 3e, pN —> eN as well 
as r —> p'y , r — 3p and the electron electric dipole moment by a numerical scan 
in both the minimal and the custodial RS model. The combination of p —> ej 
and pN —> eN currently provides the most stringent constraint on the parameter 
space of the model. A typical lower limit on the KK scale T is around 2 TeV in the 
minimal model (up to 4 TeV in the bulk Higgs case with large Yukawa couplings), 
and around 4 TeV in the custodially protected model, which corresponds to a mass 
of about 10 TeV for the first KK excitations, far beyond the lower limit from the 
non-observation of direct production at the LHC. 



1 Introduction 


Rare lepton decays are among the promising indirect probes for physics beyond the Stan¬ 
dard Model (SM). This is especially true for processes featuring lepton-flavour violation 
(LFV). For such decays the SM contribution vanishes for all practical purposes due to 
the small neutrino masses, and the experimental signature is typically very clean. The 
absence of a SM background makes the study of the decays fi —)■ ey and p —> 3e as 
well as muon-to-electron conversion in nuclei in extensions of the SM also theoretically 
cleaner than the study of lepton observables unrelated to flavour violation. The reason 
for this is that in the latter case the search for new physics involves the search for tiny 
deviations in observables such as (g — 2)^, which must then be measured and predicted 
to very high precision. 

The warped extra-dimensional Randall-Sundrum (RS) models, originally introduced 
to address the relative weakness of gravity [T] and the gauge-gravity hierarchy prob¬ 
lem 121, generically also have a rich flavour structure. Further, they provide a geometric 
interpretation for the flavour hierarchies observed in the SM: the wave functions of SM 
fermions in the fifth dimension when localized to varying degree naturally generate a 
hierarchical flavour sector j3jj6]. Despite the fact that the RS model is to some ex¬ 
tent protected from large flavour-changing neutral currents (FCNCs) |7|, the strongest 
constraint on the scale of the extra dimension comes from Kaluza-Klein (KK) gluon me¬ 
diated A F = 2 processes—notably in the kaon system [8], However, these bounds can, 
to some degree, be avoided by imposing some structure on the quark Yukawa matrices 
in the five-dimensional (5D) theory [9] or by extending the strong gauge group 10 


Nonetheless, direct searches at the LHC for KK states are no longer very promising as 
the scales that can be probed in direct production are comparably low. This is precisely 
the situation were low-energy precision observables can help considerably constraining 
the model. 

In this work we calculate the lepton-flavour violating four-fermion and electromag¬ 
netic dipole operators induced by the RS bulk fields at the tree- and one-loop level, 
respectively. We then study the p —* 3e and /j —> ej transitions, and muon-to-electron 
conversion in nuclei in the minimal and custodially protected RS model with three differ¬ 
ent implementations of the Higgs held localized on or near the TeV brane of the model. 
Lepton-flavour violating processes have been studied in the context of the RS model in 
the past, beginning with (6 11,12 


Particularly relevant are 13 , which gives the first 


comprehensive analysis of charged LFV, and 14 with the first fully five-dimensional 


treatment of loop effects, which dominate the /i —> ey decay. However, neither of the 
two can be said to provide a complete description of the loop-induced dipole operator 
coefficient, which leads us to reconsider LFV phenomenology building on our recent 
work 15 17 on penguin transitions in the RS model. In this approach we integrate out 


the five-dimensional (5D) bulk and match the RS model on the SM extended by gauge- 
invariant dimension-six operators 18,19 , which is then used to study the processes of 
interest. This two-step procedure is justified, since the scale of the fifth dimension, 
related to the mass of the lightest KK excitation, is already constrained to be much 
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larger than the electroweak scale. This approach allows for a transparent and complete 
calculation of gauge- and Higgs-boson exchange induced dipole transitions, which has 
already been applied to the flavour-diagonal case of the anomalous magnetic moment of 
the muon 15-17 . Here we focus on the flavour-changing leptonic transitions. We also 


add a 5D treatment of the bulk Higgs and the non-decoupling effect in the localization 
limit recently pointed out in 20 . Depending on the exact realization of the RS model 


loop-effects can be sizeable also for /i —» 3e and muon conversion, which at first glance 
are dominated by tree-level physics, and may introduce new correlations between /i —» ey 
and those observables. In our study we do not impose additional flavour symmetries on 
the Yukawa matrices of the model and stick to the so-called minimal and custodially 
protected models. This implies that we do not attempt here a description of neutrino 
masses and mixing in the RS framework - in the minimal RS model the neutrinos re¬ 
main massless - and only consider charged lepton flavour violation, which arises in the 
RS models independent from neutrino masses. In order to simultaneously generate hi¬ 
erarchical charged lepton masses and large neutrino mixing, RS models with additional 
flavour structure such as minimal flavour violation or discrete symmetries were consid¬ 
ered 21-23 , as well as mechanisms to suppress lepton flavour violation 24,25 . These 


have implications for LFV at tree-level, even eliminating tree-level LFV completely, but 
do not modify the structure of loop-induced LFV processes, whose complete computation 
is the main result of this work. 

The outline of this paper is as follows. The starting point of our analysis in Section[2]is 
the effective SM Lagrangian including dimension-six operators. We restrict ourselves to 
operators that can be generated in the Randall-Sundrum model at tree-level and to loop- 
induced dipole operators relevant to LFV processes. After expressing the Lagrangian in 
terms of fields in the broken phase, we match onto an effective low-energy theory. This 
allows us to use results from the literature to determine the LFV observables as functions 
of the dimension-six operator short-distance coefficients. In Section [3] we discuss three 
variants of localizing the Higgs field on or near the TeV brane and give the results for 
the corresponding short-distance coefficients in the minimal and custodially protected 
RS models using and extending results from 15,17 . LFV phenomenology of /x —> ey, 


H —> 3e and /x—e conversion and the correlations among these observables are presented in 
Section|4j We further correlate these observables with the electron electric dipole moment 
(EDM) and compare the constraining power of present limits on LFV observables and 
the electron EDM. We conclude in Section [5] 


2 From the 5D theory to low-energy observables 

In conformal coordinates the metric for the five-dimensional space-time of the Randall- 
Sundrum model is given by 


ds 2 = 



( 1 ) 
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where k = 2.44 • 10 18 GeV is of order of the Planck scale Mpp The fifth coordinate z can 
take values from z = 1/k (Planck brane) to z = 1/T (1R brane). The scale T determines 
the mass of the KK excitations in the effective 4D spectrum and is assumed to be of the 
order of a few TeV. We denote by e = T/k the small ratio of the two scales. 

The specific RS model is characterised by the particle content and the associated 
boundary conditions on the two branes. The two most common models are the mini¬ 
mal and the custodially protected model 28,29 . The latter has an extended particle 


spectrum and gauge group to enforce a protection mechanism for the electroweak p pa¬ 
rameter and the Zbb vertex. In both models the Higgs is localized on or close to the IR 
brane. The prescription how the Higgs is localized is part of the model setup itself, see 
also 30,31 . Here we will closely follow the notation and conventions of 15,17 for the 


minimal RS model and the model with custodial protection. We therefore refer to these 
references for the exact definition of the 5D Lagrangians. 

However, it is useful to introduce the parameters of the 5D Lagrangian relevant for 
the study of lepton-flavour violation. The 5D Lagrangian contains two parameter sets 
related to flavour. On the one hand there is the 5D mass My, for each 5D fermion 
held ijji. The corresponding dimensionless 5D mass parameter is c^ i = M^Jk. To 
obtain a phenomenologically viable low-energy theory that reproduces the measured 
lepton masses the cp i have to take values not too far from 1/2 for SU(2) doublets and 
from —1/2 for singlets. On the other hand, the 5D Higgs Lagrangian incorporates 5D 
Yukawa couplings t/ 5D L It is customary to work with the dimensionless Yukawa matrices 
Y = t/ 5D ) k in the minimal modelj^] and Y, Y u in the custodial model. We generally 
assume that each is roughly 0(1) and anarchic. Furthermore, the Higgs held is localized 
on or near the IR brane. There are several ways to do this. A genuine brane Higgs is 
exactly localized on the IR brane and its potential 5D structure cannot be resolved by 
construction. An IR localized bulk Higgs can be resolved, and one also needs to consider 
possible effects of Kaluza-Klein states of the Higgs as was pointed out in 20 . If the Higgs 


is exactly brane-localized one can in principle introduce a separate Yukawa coupling for 
the coupling of right-handed doublets and left-handed singlet to the Higgs. However, 
for simplicity we assume that these so-called wrong-chirality Higgs couplings 13,32 are 
equal to their “normal” counterparts. In the following, we only discuss the different 
Higgs localizations separately if they give rise to distinct results. 

Ultimately, we match the full 5D Lagrangian onto an effective 4D dimension-six La¬ 


grangian following 15 and determine the physically interesting observables as functions 


of the corresponding Wilson coefficients. The matching calculation is performed in Sec. [3] 
In the remainder of this section we take the SM effective theory including dimension-six 
operators as the starting point to compute the low-energy observables. We will con¬ 
sider only operators of dimension six that can be generated by the RS model at tree- or 
one-loop level 0 


1 For the exactly brane-localized Higgs field. For the bulk Higgs case, the relation between the 5D 
Yukawa coupling Y& and Y is given by (421. 

2 For recent model-independent analyses of LFV using the effective Lagrangian, see 33 34 . Note 
that these references use different conventions for operator normalisation, covariant derivatives and 
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2.1 The dimension-six Lagrangian 


Let us start by considering the effective SU(2) xU(l)-symmetric Lagrangian up to di¬ 
mension six, 


C = C 


SM 


JL i • /i 

£ami-o _|_ higher-dimensional operators , 
T z 


( 2 ) 


where we did not include the single dimension-five operator as it cannot contribute to 
charged LFV. We also extracted the scale, where the effective description breaks down, 
given by T in the RS model. At dimension six the operators relevant to the following 


analysis of LFV observables can be taken directly from 18,19 


£?“v = afAL^E^B^ + a^Al^a^E^W* + h.c. 


b^YL^L^L,) + bf^UYL^E^E,) + b^(E,YE i )(E 1 YE J 


+ cl j <S>'i i d“<S>(E n „E j ) + + 4$ t ib r D'‘$(L ( T' 4 7p L,) 


+ h ij (&$)L i $E j + h.c. 


£=E,L q=Q,U,D 


where L (Q) denotes doublet and E ( U , D ) singlet lepton (quark) fields. Furthermore, 

&ir A D^ = |($t i T A D^ - iD^ t a $) and with the 

covariant derivative defined as 


D„=d„- ig'^B„ - igT A W *, (4) 

and T A the SU(2) generators in the appropriate representation (e.g. T A = r A /2 with t a 
the Pauli matrices for the doublet), and Y the hypercharge. The hermitian conjugate 
in ([3]) only applies to terms in the same line. Of all the operators given in ([3]) only the 
dipole operators in the first line cannot be generated at tree-level in the RS model. The 
restricted flavour structure of some of the four-fermion operators anticipates properties 
of the RS model at the tree-level. Even though we aim to study lepton flavour violation, 
the operators with two quark and two lepton fields have to be included in C^py, since 
they contribute to muon conversion in nuclei j^] 


momentum flow compared to the present work. 

3 It is worth noting that since the 5D RS models are non-renormalizable and themselves only effective 
descriptions of some yet more fundamental dynamics, the 5D generalizations of the operators appearing 
in § could already be present as higher-dimensional operators in the RS Lagrangian. These effects 
would be suppressed by a factor (T/Auv) 2 ^ 1, where Auv is the UV cut-off of the RS model, relative 
to the effects we consider here. Since the coefficient of the dipole operator is loop-induced in the RS 
model, the underlying assumption of our treatment is that there is no tree-level dipole effect at order 
1/Auv in the UV completion of the RS model, which could compete with the loop-induced transition 
at order 1/T 2 . 
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The transition from the theory with unbroken to the one with broken electroweak 
gauge symmetry proceeds via the standard substitution rules, 


<f> ->• 


M v+h+iG \ 


Q '" Pl[ u, fa 


and 


Dfj, —y dn 


ieQA 


l 9 /j ,3 
Cw 


’W 


Q)Z, 


( Vj\ 


(5) 

U*Pl 

Ei —> VijPnlj 

vJ 


Ui -+ V^PrUj 

D, -> V*P R d, 

(6) 

■ -^(T 1 + iT 2 )W+ - 

^(T 1 - iT 2 )W ~, 

(7) 


with Zn = c w W^ — s w B IJ} and A M = c w B^-\-s w W^. Inserting (j5| and 0 into 0 gener¬ 
ates many operators, most of which cannot contribute to the processes we are interested 
in, see [l5]| for details. We are only concerned with operators that can contribute at tree- 
level (if they are generated at loop- or tree-level in the full 5D theory) or at loop-level 
(if they are generated already at tree-level). This leaves us with the Lagrangian 


/’dimfi _, /’broken 

l -LFV l ’LFV 


qA- + CxA* 

IJ ' Ji 


V M 


V2 




A A ~k 

CL-- — CL-- 

IJ Jl 


V2 




+ Pijki(£i'Y tl PRtj)(£k'Yij,PR£i) + 

- 7.‘f-V‘(d7,.PHb) - hi + 

4 C\Y 4 C\Y 

2 

+ ll^W+''‘(u i% PF j ) + h.c. 

+ 'ln — ^hi p dj) + Vn^j^ViPdj) + hx. 


+ Pw{tiYPRtj)(u k 'i ii P R u l ) + Pf£,(lil l ‘Pdj)( d kl l ,PRd{) 
+ Pijkl hidPdi)( U PlR P L U ‘) + P!$i(iiYPRti){daRPLdi) 


+ “Pdi^u^Pgu,) + p&itffPJjXd^PRd,) 

+ Wm - f>Z£mYPdj)h-ki„Pm) 

+ (A/?, + H^SmYPdi^dk^PLd,), (8) 

with P L/R = |(1 =F 75 ), and 

<4 = \ u ' aA n ,, 7.7 = £ Fuvn, <4„, 

m,n 
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<te = E P'Ub’nnVn,, % = E [U'] im U n , <? mn <1 = 2,3), 

n,m m,n 

/J$,= E (9) 

m,n,o,p m,n 

Here = c w afj — s w a ^, <G (F, V, C/“, F u , H rf , H d } are the appropriate flavour 

rotation matrices for the fermion F^\ and a similar definition applies to . 


2.2 From the effective Lagrangian to LFV observables 

The main observables for charged lepton flavour violation are radiative transitions of 
the type £\ —> f^y, lepton conversion in a nucleus, and tri-lepton decays £\ —* l 2 l£. 
These processes are usually studied in high intensity, low energy set-ups. The typical 
energy release of the process is the mass of the initial (charged) lepton, a muon or a tau. 
Starting from the effective Lagrangian at the electroweak scale discussed in the previous 
section, we construct an effective low-energy Lagrangian by integrating out the heavy 
gauge bosons and quarks and the fluctuations associated with scales above the charged 
lepton mass. For /1 —> ey, muon conversion and p —)■ 3e this low-energy theory has been 


discussed in great detail in the literature, see e.g. 27,35 , and especially 


We follow 26 and consider first the radiative decay p —> ey. The Lagrangian takes 
the formE] 


7 =A R m p £ e a ap F CTp P R £ p + A L m p £ e a ap F ap P L £ p + h.c., (10) 


where the label on the field £ denotes the lepton flavour. This Lagrangian is supposed 
to be valid at scales below the electron mass—all quantum fluctuations involving lep¬ 
tons have been integrated out and have been absorbed into the two coefficients of the 
dipole operators. In fact instead of this Lagrangian we might just as well consider the 
general U(l) em invariant vertex function for on-shell fermions (the photon momentum q 
is ingoing) 


r P (p,p') = ieQ e u p (p',s' 


rnr) + 


, ia pu q u 2 a pv q v 2 


2m,, 


W) + 


2m, 


- 75-^3 (q 


+ (gV - ) 75^4 (g 2 ) 


U e (P,s) . 


( 11 ) 


The on-shell dipole form factors of the electromagnetic muon-electron vertex are related 
to the coefficients A R and A R by 


A r = 


Qee{F 2 (0)-iF 3 (0)) 


Am?, 

F 


Ar, = 


Qee(F 2 ( 0) + iF 3 ( 0 )) 


4m?, 

F 


( 12 ) 


4 Our convention for A r , r and gi below differs from 26 , Eq. (54) by the factor —AGf/V% and 
complex conjugation. 
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where Qi = —1 is the electron charge in units of the positron charge e. Up to terms 
suppressed by powers of the electron mass the branching fraction can be written as 


m, 


Br(/we 7 ) = —+ 

47rr„ 


(13) 


Here T M is the total decay width of the muon. The generalization to l % —>■ Ifi is obvious. 


The process [i —> 3e is described by the extended Lagrangian 26 


>3e A R m/ e a ap F ap P R ^ + Atm il i e o ap F ap P L i ll 
+ gi lePdti lePde + 92 l e Pdfi 4-Pl4 
+ 93 lei'PRlp lei Me + 94 leluPrfe 

+ 95 lel V Pd» lelvPde + 95 hY Pd» ll.Pde + h.C. . (14) 


Note that the coefficients A R R and g t have mass dimension —2. The appearance of 


the same coefficients Al,r as in (10) indicates that all quantum fluctuations are again 
integrated out. In practice, absorbing e.g. electron loop diagrams involving a four-fermion 
operator into Aj j R and into a loop correction to the gi is convenient as we do not have to 


treat the different lepton flavours separately. In particular, in writing (14) the effect of 


the off-shell (q 2 ^ 0) form factors in ( JTTj ) is absorbed into the g t coefficients (see the 26 
for details). In any case, since this represents a loop correction to the Wilson coefficients 
which are already generated at the tree-level, we neglect these effects in our calculation. 

The branching fraction of fi —> 3e can easily be expressed through the coefficients g t 
and A LtR [26] p] 


Br (fj, —> 3e) = 


m, 


15367t 3 T 


v L 


\gd + \92 


+ 2 (| ^?31 2 + liul 2 ) + l^l 2 + \g& 


— 8e Re [A R (2g\ + g* & ) + A L (2g^ + g§)] 


+ 64e 2 (ln^-(l)(|yl 1 | 2 + |/l R | 2 ) 
m P 8 


(15) 


with Y p the muon decay width. The first line arises from tree-level KK exchange in 
the RS model, while the second and third involve the loop-induced dipole operator 
coefficients. The reason for keeping these formally suppressed terms is not only the 
logarithmic enhancement paired with a large numerical coefficient. One-loop corrections 
to the gi may have a similar logarithmic enhancement and the large factors of 8 and 
64 might be misleading, since it can be absorbed into the definition of the coefficients 
Al,r, which come with a 1/(47 t) 2 loop suppression. The important point is that one-loop 

5 The sign of the interference term (second line in ( fl5] ) ) depends on the convention for the covariant 
derivative. In the convention of [26] the sign is ‘+’. This is compensated by the Wilson coefficients 
Al,r, the sign of which is also convention dependent. 
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corrections to the g t introduce only a small shifts in the coefficients without altering the 
general properties. On the other hand, the dependence on A RR implies sensitivity to 
different aspects of the underlying model. In RS models A R>R have a specific dependence 
on the Yukawa couplings, which can provide an important contribution to the branching 
fraction in sizeable parts of the model parameter space. In these regions the effect of 
Al,r should not be neglected, as it will significantly alter the signatures of the RS model 
in flavour observables. 

Muon conversion in nuclei is mediated by both, operators containing quark fields and 


electromagnetic dipole operators. The effective Lagrangian is 26,36 


A R m/ e a ap F ap P R £ p + A L m/ e o ap F ap P L i p 


+ 


q=u,d 

E 

q=u,d,s 

a.m 


c q VR £el Up R^qivq 


q=u,d 


c q VL LYPd^qivq 


q 0 TD P - 


q=u,d,s 


Ml CsR ^ L 


_ML 

M 2 H 


t e P R t p G A '° p G A ap 


+ - 7 -^dc R i e P L £,,G A ’ rTp G A + h.c. 


Ml 


gg 


ap 


(16) 


Here M R denotes the Higgs mass, and G A U is the gluon field strength tensor. We do 


not include operators with pseudo-scalar, axial vector or tensor quark currents. Their 
contributions are suppressed by the nucleon number of the target nuclei and can be 
neglected. We also neglect the strange quark in the vector operators, since the coefficient 
is not enhanced by the strange-quark mass. The conversion branching fraction depends 
on properties of the nucleus that participates in the reaction. The expression, taken 


from 35 36 and adjusted to match our conventions, is 


Br (/iN —» eN) = 


m, 


4T 


capture 


A r T> + 4 


^(cL-12 *<%*,) S’+ <5^ 

H 


+ {p —>• n} 


+ {L i?} , 


(17) 


where T capture is the total muon capture rate for nucleus N. The coefficients V, V p ^ n , S p ^ n 
(the superscript refers to the proton and neutron) encode properties of the target nucleus, 
see |35l. The tilded coefficients are defined as 


rp 


^ C SLfq ) 

q=u,d,s 


(18) 


Cl, gs 


r L f p 
L ggJ Qi 


r<p 


E 

q=u,d 


r q r 

L VLJVg ■ 


(19) 

( 20 ) 



















Figure 1: Left: Matching of LFV couplings of the Z-boson onto four-fermion operators in 
the low-energy theory. Right: Higgs exchange diagram that contributes to g\ j2 and c q SR , L . 


and analogously for the p —> n and L —$■ R cases. The form factors ff’ n and ffP 
parametrize the coupling strengths of the quark scalar and vector currents of flavour q 
to nucleons, respectively, /£’" represent the scalar couplings of heavy quarks (c, b or t). 

We now determine the coefficients of the low-energy effective Lagrangians in terms 
of the Wilson coefficients of the dimension-six Lagrangian. The tree-level matching 
of the four-fermion operators in & to those in (14), (16) is straightforward. Further 


contributions arise from the lepton-flavour violating Z-boson interactions ici^PL^RfjZ^ 
together with a SM coupling once the intermediate Z is integrated out, see Figure [T| 
(left). For example, in case of fi —* 3e, the insertion of I^PrI^Z^ evaluates to 


iM = ijl 2 


2 s 2 w - 1 


u e (p 2 )'y t "PRU fl (p 1 ) u e (p 4 )'y fl P L v e (p 3 ) 


SwU e (p 2 )YPR u ^Pi) Ue{pA)lnPRVe{pz) + Fierzed diagram 


( 21 ) 


which gives a contribution to g 3 and g§. Similarly, insertions of y 2,3 will lead to con¬ 
tributions to #4 and ge. Contributions to Cy R follow analogously. Thus we fold the 
relations: 


9i — 92 — 0 ( 22 ) 

93 = ^2 ( s wl l i2 + @1211 + ® ( 23 ) 

94 = Tpi ^—^2 -(T 12 + Pi 2 ) + ft 121 1 + @\\\2^j (24) 

^5 = ^(^^7 1 1 2 + ^ 2 ) (25) 

96= ^ ( 4 ( 7i 2 2 + 7i 2 ) + Piii) , (26) 


r u — 
C VR 


1 

2T 2 


@1211 + @1211 + 2^12 



(27) 
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(28) 


C VL ~ :'2'I 1 


01211 0 1211 + ftm + 0 (712 + 712 ) (l q S 


.2 

’W 


Wr 2X" 2 

c d = — 
c 'V / l 2X" 2 


ftai + /^12H + 0 7l2 ( 1 + o S 


dEQ 


1,1 


’W 


^mi + + 2 ( 7 i 2 + 712) + 2 S w 


r 9 - 
C SL — 


V2m u T 2 


rj\2 


c q - 
C SR ~ 


s/2m u T 2 


M]l2, 


as well as 37,38 


c L = 
99 


E 


12tt — 

q=c,b,t 


C R = __ 

99 12vr 


E' 

q=c,b,t. 


-SR ' 


(29) 

(30) 

(31) 

(32) 

(33) 


It should be noted that g\ and g 2 receive contributions from the tree-level Higgs 
exchange diagram, Figure [l] (right diagram), with an insertion of one flavour-changing 
Higgs operator hliPnlj + h.c., but these are suppressed by powers of the electron mass 
(light lepton mass in the general case) and we neglect them. The same diagram (with 
the two fermion lines on the right being quarks) also generates c q SL . Here we should 
comment on a (well-known) subtlety. Naively, the operator Li&Ej&Q in (J 3 ]) modifies 
the Yukawa interaction according to 


Vij 

V2 


hZiPn&j —> 


Vij 

V2 


htjPj-f ([ j hjj 


3v 2 

2y/2 T 2 


hliPnlj 


(34) 


after electroweak symmetry breaking but before flavour rotations. Here y tJ is not the 
SM Yukawa coupling but the coefficient of the operator L^Ej in the dimension-four La- 
grangian. However, the fermion mass matrix is also modified by dimension-six operator, 


m v , „ _ vn v h y3 

s/2 s/2 3 2 s/2T>' 


(35) 


Since the flavour rotation matrices U and V by construction diagonalize the modified 
mass term m^, we have to rewrite the shift of the Yukawa couplings as (see also 1321) 


{y/2 Vij ^2 V2T*) 


M t P r £j — > 



— h 


v 


V2T 2 


hliPnlj . 


(36) 


As a consequence the factor 3/2 in the flavour-violating Higgs interaction htiPftlj + h.c. 
of ([S]) must be replaced by 1 for the computation of c q SL , c q SR above. 

We did not include the effect of the RS bulk on the effective gluon-Higgs interaction. 
The operator OLghG^G^ receives sizeable contributions from KK fermions in the loop 
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Figure 2: Example of a Barr-Zee type diagram. The box denotes the insertion of the 
lepton-flavour violating Higgs interaction. The internal gauge boson can be a Z or a 
photon. 


30 

31 

39 

40 


These corresponding contributions to c gg R after integrating out the Higgs 
loson are formally of the order u 2 /T 4 , but are enhanced by the traces of products of 
5D Yukawa matrices. We only work in leading order in the 1/T 2 expansion, so we drop 
these terms based on power counting arguments. In practice, it turns out that c g ^ R 
(with or without this additional correction) gives a much smaller contribution to the 
LFV branching fractions than, e.g., Cy R , and could be ignored altogether. 

The determination of the coefficients Al,r is more complicated. One can identify 
three contributions: (1) from tree or one-loop diagrams involving the operators in the 
dimension-six Lagrangian 15 10 - (2) from dimension-eight operators, which may be¬ 
come relevant if the dimension-six contributions are suppressed. We will discuss them 
later specifically in the context of the RS model. (This contribution can effectively be 

(3) from enhanced two-loop 
. , see 42 for a 

discussion in the context of /i —y ey. An example diagram, which avoids the coupling 
of the Higgs boson to a light lepton through the coupling to a top or gauge-boson loop, 
is shown in Figure [2} These terms are known to give sizeable contributions in models 
where the Higgs interactions are the dominant sources of new flavour violation. In the 
RS model this is generally not the case. Nonetheless we include these terms as they may 
become relevant in specific scenarios. We obtain^ 

Q t e 


included via a modification of the a R w Wilson coefficients. 

“Barr-Zee type diagrams” with a flavour-changing Higgs coupling 


m u A R = a 


12 




- E 


fc=l,2,3 


16tt 2 T 2 


m £k p 


EL 

lkk2 


3(4^) 2 r 2 I 771 ^ 7 * 2 + ^ l2 > + me7 y + m ^ 2 ~~ - \™ e i\ 2 


+ A B z 

tti^Al = [a A] }i2 


hi 2 


V2T 2 _ 

lQrt 2 T 2 "° tkHk21k 


(37) 


Q e e 


v / 2T : 




EL 


k= 1,2,3 


6 In practice, we can drop the terms proportional to the lighter lepton mass, here m e . 
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3 (ff r2 ( S W \ m e(ll2 + 7l 2 ) + m M7l 2 ] + m ell2 ~ ~ 7^ 


+ A 


BZ 


7i 2 


,2 1 


\/2T 2 


(38) 


where 42 


^4_bz — 


Qiea em \/2Gpv 
32^ 


23 3 

2 N c Q 2 J(r t ) - 3 f(r w ) - —g{r w ) - -/i(r w ) 


/(rw) “ 3(rw) + 1 2N c Q t f(r t , r„) 


2 r w 
,2 / 2 


4s 2 


4r 2 


1 ~ 1 3 

- 2 ( 5 - Sw/ c w)f(. r W i r wz) - g ( 7 _ 3sw/cw)g(r w ,r W z) ~ -^g(rw) 


-h(r w ) - -7 ^ ^ ( /(nv, ?Wz) — g(rw,rwz) 

4 4r w 


1 


4s 2 y ~ e 


Df a \r w ) + + Df c \r w ) + + Df e \r w ) 


Df a \r z ) + £>( 4 b) (r z ) + ^ 4 c) (rz) 


(39) 


with r t = m t 2 op /M 2 iggs , r tz = m 2 op /m|, r wz = m 2 w jm\ , ryy = m^/M 2 . ggs , = 
m|/M^ iggs , and iV c = 3 ,Q t — 2/3. The functions f,g,h and f,g can be found in 


42 


and the functions D* in |43|. Note again that relative signs in (37), (38) depend on the 
convention for the covariant derivative. 

All expressions for /i —> 3e and /i —> 67 in this section can be trivially extended to 
r —> ey or r —>• 3/i by exchanging the appropriate flavour indices, masses and widths. We 
note that we do not take into account the running of the Wilson coefficients between the 


high and the low scale, see e.g. 44-46 for the anomalous dimensions of the dimension-six 
operator basis. 


Wilson coefficients for lepton flavour violation in 
the RS model 


Several of the dimension-six operator Wilson coefficients have already been computed 
in 15 for the minimal RS model and in 17 for the custodially protected model. These 


works focus on the anomalous magnetic moment of the muon and therefore only included 
a subset of the operators that are needed for studies of lepton flavour violation. We 
extend the computation to the flavour-off-diagonal transitions here. 


12 























3.1 Treatment of the 5D Higgs field 

Before going into the details of the determination of the missing Wilson coefficients 
we need to address the treatment of the Higgs held. It is now well established that the 
physics of the RS model with an IR-brane localized Higgs depends on how the localization 
is implemented, see e.g. (l5][l6,30,3l] 39|. In effect one needs to specify whether the 5D 
structure of the Higgs held, that is, its 5D wave function, can be resolved within the 
model or not. If one regularizes the delta function in the fifth coordinate z, which 
localizes the Higgs close to the IR brane, by a narrow box-shaped prohle, 

S(z - 1 /T) = lim jO(z - , (40) 

5—>0 0 1 


this is equivalent to specifying the order of limits for the regulator 5 and the regulator of 
the 4D loop integrals, for example a dimensional regulator (e —> 0) or a cut-off (A —y oo). 
Removing the regulator e hrst while keeping 6 finite corresponds to the case where the 
Higgs localization width S/T can be resolved by the modes propagating in the loopj^] 
If <5 —>■ 0 hrst, the Higgs width remains unresolved. The latter case corresponds to a 
truly brane localized Higgs, whereas the hrst scenario assumes a “narrow bulk Higgs” 
(in the terminology of |3lj). Due to these non-commuting limits the RS model is only 
fully defined with a prescription on how the order of limits should be taken. 

When the Higgs held permeates the bulk, even if only close to the IR brane, the 
question arises what is the effect of the KK Higgs states. Since their masses are of or¬ 
der T/S these effects have previously been assumed to be small and have not received 
much attention. However, recent work 20 has shown that the sum of all KK Higgs 


contributions does not decouple in the localization limit. In the following we include a 
computation of Higgs KK effects in the 5D formalism. In order to have a consistent de¬ 
scription of KK Higgs modes, we abandon the ad-hoc regularization (40) and implement 
the Higgs held as a full 5D scalar doublet. Since we require both a vacuum expectation 
value and a zero-mode prohle that is strongly localized near the IR brane we need to 
introduce additional brane potentials on both the IR and UV brane. We will follow the 
construction of [47] (see also 48]) and thus use the same Higgs prohle as in 1201. The 
details of this realization together with useful formulae are collected in Appendix [Aj The 
5D prohle of the vacuum expectation value takes the form 


v(z) = yww k^T^v SM (41) 

where r?sM = v ~ 246 GeV denotes the SM Higgs vacuum expectation value (vev), and 
the zero mode prohle is, up to small corrections of order r; 2 /T 2 , proportional to v(z). 
The parameter /3 is related to the 5D mass of the Higgs held and determines the degree of 

'This is more readily seen with a momentum cut-off A. For A —> oo with fixed small but finite S 
the particles propagating in the loop can resolve the Higgs localization as the loop momentum can be 
larger than T/S. 
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IR localization; the larger (3 the stronger the localization. Since we start with a genuine 
bulk held it is always implied that f3 is finite until all other regulators have been removed. 

In order to obtain the correct SM parameters in the low-energy limit the Yukawa 
matrices and the Higgs self-coupling must themselves depend on (3. For the Yukawa 
matrices we indicate this dependence by a superscript f3, while no superscript refers to 


the /^-independent, dimensionless matrix. The relation is (see Appendix A.3) 


Y^ = 


Y 2 — CLi + CEj + (3 
\fk a/2(1 + (3) 


(42) 


with Cl., CEj the 5D mass parameters of the lepton fields in the Higgs-Yukawa interaction. 

Ultimately, we are interested in large values of (3. Whenever we give a result for the 
bulk Higgs case that does not show an explicit dependence on /3, we tacitly assume that 
the (3 —>■ oo limit has been taken, and the result should be valid up to corrections of 
0(1/(3). We will consider three different implementations of the Higgs held: 


an exactly brane localized Higgs, that is we use (40) (necessary to avoid ambiguities 
in the calculation), but take 6 —> 0 first. 


a delta-function localized narrow bulk Higgs, that is we use (40), but keep 8 finite 
until all other regulators are removed, then 8 —> 0. No Higgs KK modes are 
considered. 


a true bulk Higgs with the /3-profile (41) and KK modes. 


The second scenario is somewhat inconsistent as a Higgs held with a resolvable width 
should be accompanied by resolvable KK excitations. We will still consider it, as it 
turns out that this precisely captures the effect of the bulk Higgs zero-mode in the IR- 
localized (3 —> oo limit of the third scenario. Below we will discuss explicitly the different 
localization prescriptions only if they lead to a difference in the Wilson coefficients. 


3.2 Tree-level dimension-six operators 

3.2.1 Four-fermion operators 


The tree-level diagram contributing to the matching of the Wilson coefficients of four- 
fermion operators is shown in generic form in Figure |3j The exchanged particle could be 
an off-shell KK gauge boson or a KK Higgs excitation. The latter vanishes for [3 —> oo 
and can safely be ignored, which can be verified by explicit analytic calculation, see 
Appendix |A.4 The contribution from the remaining gauge-boson exchange diagram 


can be inferred from known results 15,17 by adjusting hypercharge and weak isospin 


factors. In case of b 1 ^ and b l j/ E , there are two contractions giving rise to an additional 
“t-channel” diagram. In the following we summarize the results for the four-fermion 
operators appearing in 
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Figure 3: Generic topology of 5D diagrams that give rise to the four-fermion operators 
upon integrating out the exchanged particle. External states can be doublets or singlets. 
Consequently the intermediate boson can be a B or Z x or the SU(2) W gauge boson, 
if all external states are doublets. Due to the chirality of the external states the fifth 
component of the boson cannot propagate. 


The Wilson coefficient of the four-lepton operator (Ljy^Lj) ( Ej^Ej ) is given by 
Y l Y e 


b LE = 

*3 


g' (6 0 + bi(c Li ) +bi(-c Ej ) + b 2 (c Li ,c Ej )) + (g 2 - g' )b 2 (c Li ,c Ej ) (43) 


with 


bn = — 


4 ln(l/e) ’ 

,, ,, 1 (5-2c)(l-2c) 

, ' l(c) = “4 (3-2c)* I 


,2c—1 


C 2c—1 


v _ 1 (1-2c l )(1 + 2 c e ){3-c l + c e ) , 1 

tb(CL, C E ) — — - —-r-—- t- 7~rz -r 111 - 


c 2C£-1 


c - 2 c E ~1 


2 (3 - 2 c l )(3 + 2 c e ){2 - c L + c E ) el - e 2 ^" 1 1 - e " 2 ^" 1 


(44) 

(45) 

(46) 


As in 15 we drop terms suppressed by the tiny ratio e — T/k. The result for the minimal 
RS model can be obtained from this expression by setting the coupling to the Z x to zero, 
which corresponds to removing the ( g 2 — g' 2 ) term in (43). Using these expressions the 
Wilson coefficient of the operator (Efl^Ej) takes the form 


-> -C E .) ■ (47) 

Here and above Y E and Y L are the hypercharges of singlet and doublet lepton held, 
respectively. For the operator (L^Lj )(L k ~/ IJi Lj) there are contributions from abelian 
Z x or B bosons as above, and additionally from the exchange of a W boson. The 
abelian contribution due to Z Xl B exchange is given by 

Y 

&Sm,b+z = 2 ^ S iAi bi k E (c Ek -> -c L J (no sum over i, k ) (48) 

The non-abclian bosons generate the operator (Lj-r^y^Lj), which is not part 

of our basis, and has to be rewritten using the SU(2) Fierz identity 

(L i r A 'fL t ) (h^Lj) = 2 (La^Lj) (Lj^Li) - {L^U) (L^Lj) . (49) 
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We then find the Wilson coefficient of (Lj 7 M Lj)(Lfc 7 M L;) to be 

h i 3 L ki = &5w,b+z + T ( b ° + + M c ^-) + b ^ CL ^ “%)) 


~ (60 + 6 i(c Li ) + 6 i(c Lfc ) + b 2 {c Li , -c L J) . 


( 50 ) 


The Wilson coefficients of the seven quark-lepton four-fermion operators are even 
simpler to compute as there are never two identical fields and all operators but one, 
(LiYr A Lj) (Qk^nfuQi), are generated via the exchange of an abelian gauge boson. 
The result is 



lL tQ 
U ij 


Y t Y q , 2 


bo + 6i(SfQj + bi(s q C qj ) + b 2 (s£C£ i , 



Y P Y X 

j 2 - (S’ - g' 2 )h(s,C,„ —SqC qj ) 


^-(^0 + bl(CQq) + &i (c Li ) + ~ c Qj)) 


(51) 

(52) 


with f G (L, E ) and q G (Q, U,D). Sf is —1 for a singlet fermion / and +1 for a doublet, 
1/ is the hypercharge of fermion /, and Y x = T|—4tan 2 ©w/(3(l—tan 2 ©m)) with T\ = 
{—1,—2,0} for q = Q,D,U. The second line in (51) is only present in the custodially 
protected model. The dependence on the 5D mass parameters of the quarks shows 
that muon conversion depends not only on the model parameters of the lepton sector. 
However, ultimately we only need operators which are built of light quarks fields after 
EWSB, and of these only the quark-flavour diagonal part. Since both the up- and the 
down-quark sector masses are hierarchical, the RS Froggatt-Nielsen mechanism generates 
hierarchical flavour rotation matrices in the quark sector (see e.g. |49]). Consequently, the 
b 2 {c x , c y ) terms—the only terms that are simultaneously sensitive to 5D quark parameters 
and contribute to the flavour-non-diagonal lepton couplings—are suppressed for light 
quarks, and we neglect them. The only unsuppressed sources of LFV are then the terms 
bi(c L J or &i (~c Ei ). 


3.2.2 Higgs-Fermion operators 


G4 


The tree-level matching coefficients of the Higgs-fcrmion operators ) 

follow from the diagrams in Figure [4j where the ones with an external gauge field are 
related to those without by gauge invariance. 

The diagrams in the first row of Figure [4] have already been computed in 


15,17 


for the minimal and custodial RS model. Their contribution to the Wilson coefficients 
c“- = c?6ij, a = 1,2,3 is given by 


c] = 


g' 2 Y E 


1- 1 - 

r(1 + 2 ceJ(5 + 2 C E .) 2(1 + 2 c E .) In 1/el 


V In 1/e 

(3 + 2c Ei ) 2 (3 + 2c e> ) J 

1 - e- 2 ^- 1 ) 
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Figure 4: Generic topologies that contribute to operators of the type 4hiD / /I> (^ 7 ^^). 
External fermion states can be either E or L. Intermediate and external gauge bosons 
can be abelian or non-abelian, the external Higgses are indicated by dashed lines. 


(9 2 - 9' 2 )Ye 


(1 + 2 c Ei ) hi 1 /e 
(3 + 2c Ei ) 


2 ce —1 


1 _ e -2c Bi -l ’ 


(53) 


Ci = 


9' 2 Yl 

1 - 1 - 

r(l-2c I ,)(5-2c t ,) 2(1 — 2c Li ) In 1/el 

e 2 ^” 1 \ 

8 

In 1 /e 

(3-2 c Li y (3-2 C L .) 

1 - e 2 ^- 1 J 


1 9 2 - 9 ,2 )Yl 


(1 - 2 clJ In 1 /e 


3-2 c 


Li 


, 2 c ,.,-1 


1 - e 2 ^- 1 ’ 


(54) 


cl = 


1 - 1 - 

[(l-2cO(5-2cO 2(1 -2ci.) In 1/e] 

e 2 ^- 1 \ 

\ In 1 /e 

(3-2 c Li y 3-2 c Li 

1 - e 2 ^- 1 J 


As in the case of the four-fermion operators the minimal RS model results can be obtained 
by removing the terms proportional ( g 2 — g' 2 ). The Wilson coefficients are independent 
of the Higgs localization provided the limit /3 —» 00 is taken in the bulk Higgs case. 

The diagrams in the second row of Figure [4] also exist, but it turns out that they 
are numerically small compared to the previous contribution. Hence, we only give the 
explicit expression for the minimal RS model: 


K = - ^g E MT) gEl (l/T)F(c Li )YlY k , 

K = K = 5 ^ h,(l/T)fL,WT)F(-c Et ) Y tt Y^ 

with 

k* (1 + 2c) + (3 - 2c)e 2_4c - (1 + 2c) (3 - 2 e)e 1 ~ 2c 
(1 + 2c)(3 — 2c)(l — e 1 -^) 2 


(56) 

(57) 

(58) 
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Figure 5: Diagram topology that gives the dominant contribution to the operator 
($t tyL&Ej. 


A similar expression is found in the custodially protected model. The smallness of this 
contribution arises from the zero-mode profiles of the light external leptons. We ignore 
the Yukawa contributions <5c“- in the subsequent analysis. 

3.2.3 Yukawa-type operators 

The dominant contribution to the Wilson coefficient of the dimension-six Yukawa-like 
operators is generated by diagrams of the type shown in Figure [ 5 ] In the 

minimal RS model there is only one diagram as the two intermediate fermions must be 
a doublet and a singlet lepton. In the custodially protected model both triplet fermions, 
T 3 and T 4 , can substitute the singlet. The contribution to the Wilson coefficient is then 
given by 

hij = Y X ^^(imlYY'Yhg^il/T) (59) 

where N cs equals one in the minimal and two in the custodially protected model. 

For completeness we remark that the diagrams in the second line of figure [| also 
contribute to the Wilson coefficient of through derivative terms that can be 

eliminated by the fermion equation of motion, such as lJ)Li = yij&Ej. In the minimal 
model we find 

Shij = ^ g El {yT)g E .{l/T)F{c Lk )yaY? k Y kj 

^hAl/T)f Ll (l/T)F(-c Et )Y ik Y ^ yij . (60) 

Due to the appearance of the small SM lepton Yukawa matrix y this contribution is 
tiny. This also holds true in the custodially protected model, and hence in the numerical 
analysis we neglect this term. However, in studies of flavour violation involving third 
generation quarks (notably top quarks) the contribution can be sizeable and must be 
included. 

3.3 Loop-induced dipole operators 

The dipole operators are generated by genuine 5D one-loop penguin diagrams. We dis¬ 
tinguish between two classes of diagrams—those with internal gauge-boson exchange 
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proportional to one Yukawa coupling Y and those with Higgs exchange, which involve 
three Yukawa couplings. A diagram such as shown in Figure [6] below counts as gauge- 
boson exchange, since it involves only a single Y. As we only need the electromagnetic 
dipole operator for the LFV processes under consideration, we reduce the number dia¬ 
grams needed for the one-loop coefficient a tl by imposing that the external gauge boson 
is a photon. In addition, we set the Higgs doublet in the operators L^a^ w EjB^ u , 
LiT A &cr flu EjW A u to its vacuum expectation value. The complete set of non-vanishing 

for the minimal RS model and in 


diagrams can be found in 15 
model. 


17 for the custodial RS 


3.3.1 Internal gauge boson exchange 

We start the discussion with the gauge-boson contribution. There are three different 
regions of 4D loop momentum that have to be distinguished. In the first region all 
propagators are zero-mode propagators. This is only possible for the diagrams that 
are already present in the minimal RS model as the additional fields in the custodially 
protected model do not have massless modes |7|. This region corresponds to a SM 
contribution and must be removed. As explained in [15] this can be done by subtracting 
the zero-mode from a single 5D gauge-boson propagator. 

When the 4D loop momentum is much smaller than the KK scale T and the loop 
contains at least one KK mode propagator, this propagator can be contracted to a point. 
The resulting diagram corresponds to a four-dimensional diagram with an insertion of 
a higher-dimensional tree-level operator. The leading dimension-six contributions are in 
one-to-one correspondence with the one-loop matrix elements of the operators considered 
in the previous subsection. 

In the last region the 4D loop momentum is of the order of the KK scale T. Only this 
region contributes to the dipole matching coefficient a A . We extract this contribution 
by expanding the propagators in the external momenta. This automatically removes the 
low-momentum region and prevents double-counting of the insertions of tree operators 
15 . As the amplitude of the dipole operator is proportional to we only have to 

expand to linear order in the external momenta. The remaining calculation requires 
the numerical evaluation of integrals over the modulus of the 4D loop momentum and 
several bulk coordinate integrals in the interval [l/k,l/T\. The complete numerical 
calculation was performed for the minimal RS model in |l5] and for the custodially 
protected model in [l7|. Here we use the results of [l7| as the routines used there for the 
numerical evaluation have been considerably improved compared to 15 . We perform 
the calculation in 5D R^ gauge and use the spurious dependence of the numerical result 
on the gauge-fixing parameter as an additional estimate of the numerical uncertainty. 
From this we conclude that the numerical accuracy is below 0.1%, which translates to an 
error of about 5% for the flavour-non-diagonal LFV terms after rotation to the standard 
mass basis. This is sufficient for our purposes. 

For completeness we remark that some diagrams induce a scheme-dependence of the 
dipole coefficient via finite but IR-sensitive £ x 1/e terms. This scheme dependence is 
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Figure 6: Example diagram with internal KK Higgs modes. 


cancelled by the scheme dependence of the 4D one-loop penguin diagrams with insertions 
°f Pij E i lij p|. 

In principle there could be an additional momentum region. The width of the Higgs 
localization introduces the new scales T/S and (3T for a delta-localized and a bulk Higgs, 


respectively. In 16 it was shown that the effect of this additional momentum scale 


does not contribute to the gauge-boson exchange diagrams for <5 —» 0 (or, equivalently, 
(3 —> oo) when only the Higgs zero-mode is considered. For the bulk Higgs case it still 
needs to be shown that the contribution of the infinite tower of Higgs KK modes also 
vanishes for (3 —> oo. To this end let us examine the diagram shown in Figure [6} Up to 
a constant prefactor it is given by 


d d l l 


2 r^r 


dzdxdy 


Y$S% W/i? W* 101 (!/) A| MS (;, z, y)F+ ((, z, *) Y A | M S((, a , * 


( 2n) d J i k 13 z 5 x 5 y 3 


dP 


(61) 


For the explicit expressions for the zero-mode subtracted gauge boson propagator A™ s 

The Higgs propagator A $(l,z,y), its 
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and the fermion propagator F we refer to 
zero-mode subtracted version A| MS (/, z, y), the Higgs zero-mode <&(°)(y), and the Yukawa 
coupling Y :i are discussed in Appendix [a] 

We now show that the KK Higgs contribution is 0(1/(3) and therefore can be ne¬ 
glected for large f3. The Yukawa matrix Y 13 and zero-mode profile <fb°)(r/) both scale as 
yf]3. Since the zero-mode profile is localized near the IR brane, the associated 5D coor¬ 
dinate integral over y is effectively restricted to the interval [(1 — 1/(3)1/T, 1/T] of length 
1 /((3T). Hence the y integration introduces a factor oil/ (3. The integration over y then 
compensates the factor (3 from the product Y' s $(°)(y) independent of the magnitude of 
the 4D loop momentum /. For l <C T and l ~ T, the Higgs propagator scales as 1/(3 and, 
after a change of integration variables from (x, y, z} to {y, y — z, z — x}, one finds that 
the integrand is dominated by the region where the distance z — y is of the order 1/(3 (see 


also Appendix A.4). Putting all factors together, we conclude that the integrand scales 
as 1 /(3 2 for small loop momenta, and hence the integral over these momentum regions 
also vanishes for (3 —y oo. For loop-momenta l of order (3T, we can expand the fermion 
and boson propagator for large momenta, in which case they become simple and their 
dependence on the loop momentum can readily be extracted. The Higgs propagator 
is more complicated, but it can only depend on the scale (3T and therefore scales as 
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Figure 7: Integrand as a function of the loop-momentum l for (3 — 10 (black, dotted), 
f3 = 20 (blue, dashed) and f3 = 40 (green, solid). For clarity all curves have been rescaled 
relative to the maximum of the integrand for (3 = 10. For loop momenta in excess of (3T 
the integrands show a universal l/l 2 behaviour. 


l/(/3T). We find that the product of all three propagators together with the derivative 
d/dl 2 , which counts as l/(/3T) 2 , compensates the factor l 5 ~ (/ 3T) 5 from d 4 /Z 2 ~ dll 5 . 
We are left with the two integrals over y — z and z — x. For l ~ (3T the integrand is expo¬ 
nentially suppressed for \z — x\ > l/l and \y — z\ > l/l , and hence each of the coordinate 
difference integration regions is effectively restricted to size l/{/3T). We then find that 
the total scaling of the integrand in this momentum region is oc 1//3 2 . The integral over 
dl can only compensate one inverse power of /3 and we conclude that the integral over 
the region l ~ j3T vanishes as well for f3 —> oc. For very large loop momentum l f3T 
we can expand all propagators. Now all bulk coordinate differences are constrained to 
be within about l/l (l is now the largest scale) and the 5D Higgs propagator scales as 
l/l. This ensures the convergence of the integral as the integrand vanishes as l/l 2 for 
l —> oo. The parameter j3 only enters through the integral over y , which is cancelled 
by Higgs profile and Yukawa coupling, hence the integrand is independent of j3. This 
universal behaviour allows for a straightforward determination of the contribution of the 
region l (3T: 


dl _ 1 

] fT p~pf' 


(62) 


Hence the integral over this region vanishes in the large f3 limit. Since this holds in all 
regions, we conclude that KK Higgs contribution vanishes as 1//3. 

This can be verified numerically as shown in Figure [7j The three curves correspond 
to different values of /3 (10, 20 and 40, respectively). For better visibility all curves are 
normalized to the maximum of the f3 = 10 curve. The maximum of the integrand is close 
to l ~ f3T and exemplifies the 1//3 2 scaling of the integrand in that region. For large 
modulus of the (euclidean) loop momentum the three curves lie on top of each other 
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Figure 8: Contour plot of Aij in the custodially protected RS model normalized to 
its value for cl = \ce\ = 0.7 as a function of the 5D mass parameters cl and ce for 
T = 1 TeV, k = 2.44 • 10 18 GeV. 

consistent with the /3 independent asymptotic expression. Consequently, the integral 
over l as well as the contribution to the dipole operator coefficient vanishes for f3 —> oo. 

To understand the numerical size of the Wilson coefficient a?- generated by the gauge- 
boson contribution it is convenient to factorize all terms that combine to the 4D Yukawa 
matrix before rotation to the mass basis: 

4 = = y»As • (63) 

The remaining short-distance function Aij depends only on the 5D bulk masses of the 
external fermion fields with flavours i, j and the RS scales k and T. Aij can be interpreted 
as a measure of the misalignment between the mass matrix of the lepton sector and 
the dipole coefficient a?- before rotation to the mass basis. If the Aij were all equal, 
no LFV would be generated by the gauge-boson exchange diagrams. Figure [8] shows 
the result of the numerical computation of Aij for the custodially protected model at 
the KK scale T = 1 TeV. There is a small asymmetry in the dependence of Aij on 
the bulk mass parameters of the external lepton fields, which arises from 5D diagrams 
with non-abelian gauge bosons as the W bosons do not couple equally to singlet and 
doublet fields. To reproduce the 4D lepton mass matrix the bulk mass parameter cl 
of the doublet muon (electron) has to be around 0.57 (0.66) and the masses of the 
corresponding singlets around —0.57 (—0.66), if the SM mass hierarchy is carried by 
both singlet and doublets. As illustrated in the figure the variation of Aij in this region 
is around ±(2-3) %. In an extreme case where e.g. all singlets are “delocalized” with 
bulk mass parameter ce = —0.5, the bulk mass of the doublet muon (electron) has to 
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be around 0.64 (0.8), and the variation is less pronounced. For the minimal RS model 
the dependence of A l3 on the bulk mass parameters is slightly smaller in the region 
of mass parameters relevant to muons and electrons than in the custodially protected 
model 


15 
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It follows that the gauge-boson exchange contribution a®, to the dipole 


coefficient has smaller off-diagonal elements by a factor 30 to 50 compared to the flavour- 
conserving diagonal entrie^]— the RS model has a built-in protection from large gauge- 
boson induced LFV transitions. It is interesting to note that the variation of A l3 increases 
for decreasing absolute value of both bulk masses. Since typically the absolute values 
of the 5D bulk masses decrease with decreasing magnitude of the 5D Yukawa couplings, 
a smaller absolute value of the 5D Yukawa couplings leads to more pronounced LFV 
transitions from internal gauge-boson exchange. 

3.3.2 Internal Higgs exchange 

Unlike gauge-boson exchange the internal Higgs-exchange contribution depends strongly 
on the Higgs localization. 


Delta-function localized Higgs 


We first consider the delta-function localized Higgs (40). For the minimal RS model the 


result for a delta-localized bulk Higgs without KK modes was determined in 15 to be 

Q,r T 3 


< = 


r J~'o 

192t r 2 k 4 Fe ' ) 


Q,e T 3 
192tt 2 k 4 


f^(l/T)[YYA%g ( E °hl/T), 


(64) 


where 


Fe = fSWT)Y it F(-c Et )Y^(l/TrY hj g'“'(l/T), 


t AO), 


M, 


Fl = lf}[VT)Y a gm\IT^‘Y^FMY^'aiT). 


A) 


2v"t 


,( 0 ), 


(65) 


Here the three times repeated indices k and h are summed over only once and the 
function F(c) is defined in (58). The expression in the first line arises via “off-shell 
contributions” (see 15,16|) and is suppressed by fermion zero-mode factors. The last, 


numerically dominant line (64) is not present for the exactly brane-localized setup. As 


a consequence the Higgs-exchange contribution, which involves three Yukawa matrices, 
is suppressed in the minimal model with an exactly brane-localized Higgs boson. 

The corresponding result for the custodially protected model was partially determined 
However, due to a misplaced SU(2) index the two diagram topologies (absent 
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®This factor is responsible for the larger numerical error for the flavour-violating transitions men¬ 
tioned above. 
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Figure 9: Additional Higgs-exchange diagrams in the custodially protected model. 


in the minimal model) shown in Figure [ 9 ] were missed]^] These diagrams have a non¬ 
trivial dependence on the Higgs localization. In the following we compute these missing 
diagrams. Note that the exchanged Higgs refers to the zero mode, since for now we have 
adopted the theta-function regularized delta-function Higgs profile. 

To illustrate the computation we consider explicitly the sum of the two right-most 
diagrams in Figure [9] where the photon couples to the fermion line. Up to prefactors that 
depend on the U(l) and SU(2) charges of the fermions, and after some simplifications 
the contribution to the dipole operator structure is given by the integral 


Inew ~ <53 k 12 


T 3 T 9 f A2 dl 2 
0 


-l/T 


167T 2 


'(1 ~S)/T 


dx dy dz f ( L 0) (x)Y ik Y^ h Y hj g^ (z) 


1 

x - 
2 


di 2 (l A di 2 (F Xk (l,x,y)F+ h (l,y,z))) + fd^dr^d F% k {l,x,y)d F$ h (l,y,z)) 


, ( 66 ) 


where we used the notation of 15 for the different fermion propagator functions and 
applied a momentum cut-off A to the loop integral. The fermion A" can be part of the 
SU(2)i singlets £ 2 , T 3 , T 4 while Y is part of the custodial bi-doublet £ 1 . Both propagator 
functions d~F + in the second line obey Dirichlct boundary conditions at y = l/T, 
while the propagator functions F + and F~ have Neumann boundary conditions on the 
IR brane. That is, the second term in square brackets contains wrong-chirality Higgs 


couplings 13 , whereas the first features right-chirality Higgs couplings (correct-chirality 
in the language of |20|). 

Consider first the first term the square brackets in (66). The l 2 integral can be 


carried out trivially, since it is a total derivative. Only the upper limit gives a non-zero 
contribution. Then we use that for l 2 ~ A 2 T 2 and x,y close to l/T, the fermion 
propagators can be simplified, 


F x (l, x, y) « Fy(l, x, y) oc cosh(/(l/T — x))e l< ^ y 1 ^ T - ) 0(x — y) + {x <B- y} , (67) 

which allows all coordinate integrals to be evaluated analytically. The result is quite 
involved and depends crucially on the expression e _ ^ A<s ^ T , such that it vanishes for A —> 

9 Note that for external quarks, i.e. for the calculation of the electromagnetic dipole coefficient of 
quarks, the diagrams are present even in the minimal RS model. 
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oo at fixed, finite 5, but is equal to 1 for 6 —> 0 with fixed A. 1U Hence, for a narrow bulk 


Higgs (small but finite 5, A —)■ oo first) the first term in square brackets in (66) does not 
contribute to the dipole Wilson coefficient. This reproduces known results 20,50 for 


the right-chirality Higgs couplings. On the other hand, for the exactly brane-localized 
Higgs, the dipole coefficient receives a finite unsuppressed contribution. This effect is 
only present in the custodially protected model, since in the minimal model the diagrams 


do not exist. The contribution of the second term in the square brackets in (66) can be 


computed following the approach of 16 . Here situation is exactly opposite: the narrow 


bulk Higgs leads to a finite contribution, whereas the integrals vanish when the Higgs 
width 5 is taken to zero before the regulator of the loop integral is removed, i.e. for the 
exactly brane-localized Higgs. The total result is 


T 3 


■ * T jf { S(im\YY'Y} ij gVH\IT) “exactly localized” 


,(0) 


T.i r/r < 


32tt 2 T 2 k 4 

T 3 


( 68 ) 


■ t nf { S(llT)\YY'Y} lj gVH\IT) “narrow bulk” 


,(0) 


967t 2 T 2 k 4 


The contribution of the left diagram in Figure [9] can be obtained analogously. 

Putting all diagrams together and including the SU(2) and hypercharge factors we 
find in the custodially protected RS model for the narrow bulk Higgs 


i H 


Q,e T 3 


192t r 2 k 4 

and for the exactly brane-localized Higgs 

T 3 


■ 2 ' fLhl/T)[YY'Y} v gf (1/T), 


H 


13 192tt 2 k 4 


■ (- 3 ) • f'S(l/T)[YY'Y + Y n Y<Y] ijg ^(l/T) 


(69) 


(70) 


In both equations we neglected terms similar to the first line in (64), which are suppressed 


by lepton masses and/or lepton zero-mode profiles. These terms are always subleading 
in the custodially protected model irrespective of the Higgs localization, and we ignore 
them in the further analysis. Note that afj has opposite signs for the narrow bulk and 
exactly brane-localized Higgs. 


Bulk Higgs with a (3 profile 


Next we consider the case of a bulk Higgs with the fd profile (41). The dominant con¬ 


tributions were studied in some detail in 20 and numerical estimates were obtained by 


10 In dimensional regularisation the calculation is more tedious. The second line of (66) can be written 


as a total derivative plus an evanescent term oc (d— 4). The non-commuting limits manifest themselves 
in form of factors of S e or in combinations of incomplete Gamma functions like T(—1 + e,6). The final 
result coincides with the one obtained above with the cut-off regulator. 
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Figure 10: Example diagram with an internal Higgs exchange and three Yukawa cou¬ 
plings. The internal Higgs can represent the zero-mode or the tower of KK modes. 


summing a large number of KK modesj^] Using 5D propagators the effect of the Higgs 
zero mode can be computed analytically for large /3. To see this let us focus on the 
simplest diagram, shown in Figure [TOj Other contributions can be obtained analogously, 
but may require appropriate expansions of the fermion propagators for a fully analytic 
result. For light external fermions the dominant contribution can be written as 


Zh = Q u eY? k [Y%Y\ 


t ^/3 

kj 


_dH_ [T dz dy dx dw (0), , (0 ), w)d - F+(v x w) 

(2vr) 4 Ji k™x 5 y*z 5 w 5 1)9e > [ 


xA$(l,i,z)<l l0) (w) [d F+ h (p' -l,z,y)F+ h (p -l,y,x)f\f-t) 


+ F e h (p' ~ l ,z,y)d F+ h (p-l,y,x)(f - t)Y] 


(71) 


where we chose p,p' for the incoming and outgoing fermion momentum, respectively. 
The integral over the w coordinate can be taken right away as we can set p to zero in 
the external fermion propagator: 


£(*, p, c L , c E ) = J T Yu 9 f(w) d-Ft(p = 0, x, w) *'“>(«-) (72) 

J k 

. I 1 + 2 c E / 2(1 + /3) 1 (Tx) 2+Cl €~ 5 / 2 

” VI - e 1+2cfi V 1 - e 2+2/3 2 - c L + c E + P 1 - e 

l-2c L / e 2c L -l 


11 At this point it is worth recalling that the dipole transitions are not sensitive to the UV cut¬ 
off Ayv °f the RS model. That is, whether the summation of KK modes is extended to infinity or 
truncated when the KK mass reaches Auv, results in differences suppressed by inverse powers of Auvj 
which can be absorbed into the coefficient functions of higher-dimensional operators. Equivalently, in 
the 5D treatment the integrations can be done over the entire loop momentum and bulk space without 
imposing an explicit cut-off. This is completely analogous to the standard practice of performing the 
loop integral over all four-momenta for UV-finite observables in the Standard Model, even though the 
Standard Model is most likely only valid up to a certain scale. When this logic is applied to summing 
KK modes of a narrow bulk Higgs in the RS model, since the non-decoupling effects is related to the 
new scale /3T, it is, of course, implied that T <C /3T <C Auv- Otherwise the issue of Higgs localization 
could not be separated from the UV completion of the RS model. 
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After expanding the remaining integrand for small p, p' we perform the integral over the 
photon vertex bulk position y using the completeness and orthogonality relations. We 
then find 


J„ = Q„eY? h [Y^Y* 


d^l l T dz dx „(o) 


(2t) 


10 ^5 


f j( z )£( x i Pi C L k i C E X, z ) 


X 


-l 2 d^d F+ h (l,z,x) 


(p M + p ) ■ 


(73) 


This leaves us with only three integrals over x, z and the loop momentum. 

Let us first consider the Higgs zero-mode contribution by substituting A$(l,x,z) —> 
i/l 2 x Since (3 is large but finite until all integrals have been carried out 

and all regulators removed, we can perform the momentum integral directly in d — 4 
dimensions. To this end, we switch temporarily to the mode picture for the fermion 
propagator, evaluate the integral 


/ 


dH 

(2vr) 4 (l2 



i 1 
2(47t) 2 m 2 ’ 


(74) 


and resum the mode expansion back into 5D propagators, which results in 


if = -u 


MIX P 
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2(4tt) 


-d Ft(0,z,x)(p“+p"‘) 


(75) 


Since the zero-momentum limit of the fermion propagator has a simple form, the two 
remaining integrals are elementary. The final analytic expression is lengthy and valid 
for any positive value of (3. We refrain from giving the explicit expression. However, the 


limit (3 —> oo is straightforward. After using (149) to relate the Yukawa matrices for the 


bulk Higgs to the couplings for the delta-regularized Higgs we recover the same answer 


as m 


15 


r (0) 

-H 
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iQ^e 


T 3 


96tt 2 T 2 k A 


f?KVT)Y ih Yt k Y kjg fniT) (jY + p'") 


(76) 


This observation is general: the Higgs zero-mode contribution of the bulk Higgs in the 
f3 —* oo limit is always equal to the one of the theta-function regularized brane Higgs, 
see (64) and (69). In other words, the localization limit of the bulk Higgs is independent 


of the bulk profile at finite Higgs localization width. 

We still have to determine the contribution from the tower of KK Higgs excitations 
in the loop. To illustrate the computation in the 5D framework, we consider again 
the diagram shown in Figure [lOj T he KK contribution is obtained by the replacement 
A<j,(Z, a;, z) —* A| MS (Z, x, z) in (|73[). An analytical evaluation seems difficult even for 
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/3 » 1. Before turning to the numerical calculation we shall first show that the KK 
contribution does not go to zero for large (3 despite the fact that the lowest KK masses 
are of order (3T. This confirms the non-decoupling effect found in 20 , now in the 5D 
framework. To this end we look at the different loop-momentum regions separately. 
There are two relevant scales, the KK scale T and the ffiggs localization scale (3T. This 
leads to several momentum regions that allow for various expansions of the propagators. 
The expanded forms can then either be integrated directly or at least their (3 scaling can 
be determined. 


For small loop momenta I C T we can expand both the fermion and the ffiggs 
propagator around l = 0. We can then analytically integrate the x and z coordi¬ 
nates as in ( [75| . In this region the scaling with (3 must be the same as the scaling 
of the Wilson coefficient of the four-fermion operator discussed in Section 3.2.1| and 
in Appendix |A.4| That is, for large (3 the integrand scales as 1/(3. Hence, the total 
contribution from this region vanishes for (3 —>■ oo. 


• The second region is l ~ T. For the Higgs propagator we can use the same 
expansion for small euclidean momenta as for l <C T but the fermion propagator can 
no longer be expanded. Nonetheless, d~F^{l, x, z ) does not introduce an additional 
/3 dependence in this momentum region. We recover the overall scaling oc 1/(3 for 
fixed values of l just as for l <C T. The only difference to the region with / T is 
the scaling of the integrand with the loop momentum l, which no longer is a simple 
power law. However, for l ~ T, the scaling of the integral with (3 is the same as the 
integrand, that is 1/(3, and hence the contribution from this region also vanishes 
for (3 —> oo. 


For loop momentum l of the order (3T we can make use of an expansion of modified 
Bessel functions of the form Zp((3x) and 1Cp((3x) for large (3, given by 






2 i t (3 (1 + x 2 ) 1 / 4 


g(x ), Kpifix) 



,-Pf 0*0 


2(3 (1 + x 2 ) 1 / 4 


g(x). (77) 


The exact expressions for the functions /, g and g can be found in 151,52 . Here we 


only need that /, g and g depend on (3 only via terms that vanish at least as fast 
as 1//3 for (3 —> oo and that /(x) is a strictly monotonically increasing function 
of x. Using these expansions one can show that the Higgs propagator retains the 
same 1/(3 scaling as in first two regions. Taking into account the behaviour of the 
fermion propagator for l T we fold that d A l l 2 d/ 2 d~F + (l, z, x) counts as a factor 
of dl l or equivalently dl ■ (/ 3T ). This cancels the 1/(3 from the Higgs propagator and 
leaves us with the coordinate integrals. Their counting is easier to determine when 
the integral over w has not yet been carried out. The integral over w then cancels 
the \//3 factors from the Higgs zero-mode profile and one Yukawa coupling. Every 
integral over a coordinate difference counts as 1/(3 (compare the discussion of KK 
effects in the gauge contribution). Including the two remaining Yukawa couplings, 












we find that the integrand scales as 1//3 in the region / ~ f3T. Hence the integral 
over the domain l ~ f3T takes a constant value for f3 —> oo. 


• Finally, for l /3T we expand the Higgs propagator for large momenta, since it 
is now dominated by the scale l and no longer by (3T. Consequently, the Higgs 
propagator scales as 1//, and the distance \x — z\ is limited to be of order 1/7. This 
effectively trades two powers of l/(/3T) for two powers of 1 /l compared to result 
in the l ~ (3T region, resulting in the scaling oc (3/1 2 of the integrand. The final 
integral over the modulus of / is therefore convergent and since 



1 

T ’ 


(78) 


the high-momentum region also gives a finite /3-independent contribution to the 
dipole operator coefficient. 


Since in every region the integral over l either vanishes (l <C T, l ~ T) or converges 
(/ (3T and l ~ f3T) to a constant, the contribution to the dipole Wilson coefficient due 

to the Feynman diagram in Figure [10] tends to constant for large /3 as announced. For 
large values of /3 the integral is further dominated by the high-momentum regions and 
therefore the 5D masses of the fermions enter predominantly via the external zero-modes. 

The left panel of Figure [Tl] shows the numerical result for the integrand as a function 
of the loop momentum l and demonstrates the expected inversion of the order of the 
curves for different /3 values from the intermediate to the high-momentum regions J^] 
The right panel shows the KK Higgs contribution as a function of /3 normalized to the 
zero-mode contribution in the /3 —>■ oo limit. The plot illustrates the approach of KK 
contribution to a constant. The relatively fast convergence with increasing /3 is a feature 
of the simple diagram topology under consideration. The plot shows that the KK Higgs 
contribution while somewhat smaller than the corresponding zero-mode contribution is 
of the same order of magnitude 


20 


A similar scaling analysis can be applied to all other diagrams involving KK Higgs 
modes. We will not discuss them in detail, as we have to resort to a numerical evaluation. 
In Appendix[B]we give the numerical ratio of the KK tower to the zero-mode contribution 
for each diagram topology. The final result for dipole operator coefficient generated by 
the exchange of the KK Higgses is 


<•“ = ' - 4 kk ■ /®(1/T)[yyty] ij 9 g(i/T) (79) 

in the minimal model and 

“«' KK = ifti? ' A?(i/r) [a&yy'y + eg K uyJy]y< ) (i/r) (SO) 

12 Note that the solid curve for /3 = 160 does not reach the asymptotic region of very large loop 
momentum l /3T, while /3 = 10 is on the small side for the /3 1 scaling to hold. When taking the 

coordinate integrals analytically (possible in some of the momentum regions) we encounter ratios of T 
functions such as r(6 + /3)/r(7 + /3), which scale as 1//3 for large f3 , but /3 ~ 10 is not quite large enough 
to make this manifest. 
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Figure 11: Left: Absolute value of the integrand for the diagram of Figure [TO] with 
zero-mode subtracted Higgs propagator as a function of the loop momentum l. The 
curves correspond to f3 = 10 black (dotted), 40 blue (dashed) and 160 green (solid). The 
KK scale was set to T = 1 TeV. Right: KK Higgs contribution to the dipole operator 
as a function of j3 normalized to the (3 —>■ oo limit of the Higgs zero-mode contribution. 


in the custodially protected model, which should be compared to the second line of (64) 
and (69), respectively. Here we again dropped the suppressed “off-shell terms” similar 


to those in the first line of (64). The numerical values of the coefficients are 


Akk = 0.46(0.04) 


■A-kk = 1-4(0.2) 


Brk = 0.1(0.05), (81) 


where the number in parenthesis shows the estimated error due to the extrapolation to 
(3 = oo. The sizeable relative uncertainty in comes from large cancellations among 
the various contributions to the coefficient. In the minimal (custodial) model the KK 
contribution is about 50% (75%) of the zero-mode contribution. 

Irrespective of the Higgs localization, the dipole coefficient generated by Higgs 
exchange is in general misaligned relative to the mass matrix. For the bulk Higgs case 
the numerically dominant terms scale as YY'Y in both the minimal and custodially 
protected RS model (which includes Y U Y^Y). After rotation to the mass basis this 
potentially generates large LFV transitions. For the same reason, even after the rotation 
to the mass basis, unlike the gauge boson contribution, the Higgs contribution depends 
strongly on the values of the 5D bulk mass parameters and the 5D Yukawa matrices. It 
usually increases with the magnitude of the Yukawa matrix entries. In the minimal RS 
model with an exactly brane-localized Higgs, however, a H is much smaller than for the 


bulk Higgs case (recall that in this scenario only the first line of (64) is present). 


3.4 Dimension-eight operators 

The effects of dimension-eight operators are suppressed relative to the dimension-six 
ones by a factor of 0{y 2 /T 2 ) and therefore negligible. However, for LFV observables 
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this counting can be numerically upset, as noted in [14] , since the leading dimension-6 
contribution to the dipole operator from gauge-boson exchange is suppressed by a factor 
of 30-50 due to the near-alignment discussed above and in 15,17 . Relevant dimension- 


eight effects can arise directly from dimension-eight operators and indirectly from v 2 /T 2 
corrections to the field rotation to the mass basis. 

The first class corresponds to the descendant (Ljcffi" E i )$X l <f> (A" = B, W) of the 
dimension-six dipole operator (Licr ^which after EWSB give rise to the same 
dipole vertex structure. However, the dimension-eight operator has a coefficient function 
proportional to YY^Y even for the internal gauge-boson exchange contribution, and does 
not suffer from the alignment suppression of terms proportional to Y. Depending on the 
value of T, the dimension-eight contribution may then be the dominant source of flavour 
violation. This is relevant only for the case of an exactly brane localized Higgs in the 
minimal RS model, where the contributions to the dimension-six dipole Wilson coefficient 
cubic in the Yukawa coupling due to Higgs exchange are also suppressed (see previous 
subsection). 

The second class of dimension-eight effects arises from the fact that the tree-level 
relation 
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(82) 


that defines the rotations U, V to the mass basis |3| receives corrections due to multiple 
Higgs vev insertions J^1 The diagonalization condition has the form 
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-V kn = diag {m e , m M , m T } , (83) 


cf. ( |35| . The modified U and V field rotation matrices applied to the Lagrangian (|3]) 
generate an additional source of LFV which formally enters at the same level in the 
v/T counting as dimensions-eight operators, which can be taken into account by the 
substitution 


gR cR + 
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6 T 2 


a: 


ij\Y^YYlY 


(84) 


The direct effect of the dimension-eight operators is more difficult to estimate. We 
have to evaluate the contributions to the dipole-like operators that appear at the dimen¬ 
sion eight level, i.e., 
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w - d ™- 8 (L,t a Ej) IT^/h 1 4>. 


(85) 


The computation of the electromagnetic dipole coefficient afj m 8 = cos ©iy a?’ dim 8 — 
sin ©w a7’ dim ~ 8 would require the computation of roughly 150 different diagrams in the 
5D theory for the minimal RS model alone. 

13 The square root factors arise from the explicit expressions for the lepton zero-mode profiles. 
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Figure 12: Example of a diagram contributing to the matching onto the dimension- 
six dipole operator and several related diagrams that contribute to the dimension-eight 
operator. 


Fortunately, only some of these diagrams actually contribute. For the following we 
consider only the minimal RS model with an exactly brane-localized Higgs. For the 
other Higgs localizations the dimension-six dipole is always dominant and dimension- 
eight terms are negligible as discussed above. We then have two fundamentally different 
contributions: from the so-called wrong-chirality Higgs couplings (WCHC) and from the 
ordinary Higgs couplings to lepton modes with the same chirality as the SM zero modes. 
It turns out that for the exactly brane-localized Higgs the WCHC contribution can be 
computed analytically and is simply given by 


dim-8,WCHC _ g I 

a ij ~ 3 a ij ly-4-yyty 


( 86 ) 


in terms of the dimension-six gauge-boson exchange contribution. 

To illustrate how this result arises let us consider the left-most diagram in Figure [12] 
(W8 in the notation of 15 ), which contributes to the matching of the aJJ coefficient. 
There are 10 ways to add two additional external Higgs lines to the fermion line. However, 
since 8/T (5 being the Higgs localization regulator) is much smaller than the dimensional 
regulator or, equivalently, than the inverse loop momentum cut-off, we find that only the 
three diagrams shown to the right in Figure [12] give a non-vanishing WCHC contribution 
for <5 —>■ 0. In each case the integrals over the Higgs vertices can be taken analytically. 
In the above example the WCHC contributions of the two right-most diagrams cancel, 
and the remaining diagram can be expressed in terms of the associated dimension-six 
diagram as shown in (86). Similarly the descendants of all other dimension-six diagrams 
can be shown to satisfy (86). 

Hence the effect of the WCHC can be included via the redefinition 


a ij a ij ~ 


6 T 2 


a: 


ij\Y- s-yrty 


(87) 


where we used that the Higgs fields will assume their vacuum expectation value —> 
v 2 / 2 ). Combining this with (84), we find that the direct and indirect contribution 
cancel. That is, at the dimension-eight level the WCHCs do not generate sizeable flavour¬ 
changing transitions by lifting the misalignment suppression and can be ignored. 

This leaves us with the dimension-eight contributions that have no WCHCs. In the 
minimal model as defined in 15 there are no such contributions from the diagrams with 
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non-abelian vertices. Then there are only seven non-vanishing diagrams that involve 
an internal W boson, but about 50 diagrams with a hypercharge boson. Fortunately, 
the limited particle content of the minimal model allows us to recast the expressions of 
all diagrams in the form of the original dimension-six diagram with modified fermion 
lines. For instance, the second diagram in Figure 12 has terms without WCHCs, but 
differs from the original diagram only by the two additional (zero-momentum) Higgs 
insertions that modify one fermion propagator. This can easily be calculated as the 
Higgs vertices can be treated analytically. Since the flavour-dependence of the fermion 
propagators (excluding zero-modes) is relatively mild, one can use the single-flavour 
approximation, where the Yukawa matrices are the only flavour-dependent quantities. It 
is then straightforward to compute the contribution to the dimension-eight coefficients. 
We find 


dim—8 


—n a I 

U u ij |y_>yyty 


( 88 ) 


This size is in agreement with the estimate given on the basis of a subset of diagrams in 


14 


where the non-abelian contribution was found to be 8 
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The minimal model requires a KK scale T > 4 TeV in order to pass the constraints 
set by electroweak precision observables [49). Hence the dimension eight contribution to 
ctjj is suppressed by an additional factor v 1 /T' 2 of at least 1/500. We therefore neglect 
the contribution of the dimension-eight terms to the off-diagonal elements of ay,-, since 
it is smaller than the effect of the Barr-Zee diagrams which also feature three Yukawa 
couplings without the need to take the dimension-eight term into account. 

For the custodially protected model the dimension-eight coefficient would be much 
harder to compute. Not only does the number of non-trivial Feynman diagram topologies 
increase significantly, but the larger particle content leads to numerous non-vanishing 
possibilities to assign the various fermion species to each topology. However, independent 
of the Higgs localization there always exists an unsuppressed dimension-six contribution 
proportional to YY'Y, hence the dimension-eight terms are never relevant. 


4 Phenomenology 


The Standard Model in its original form 53 does not allow for flavour violation in the 


lepton sector. Even after introducing neutrino masses, the charged LFV processes are 
suppressed by the tiny neutrino masses and too small to be detected in any foreseeable 
experiment. Any signal of charged LFV is a clear sign of physics beyond the SM. In our 
analysis of LFV in RS models, we focus on the three processes with the highest current 
experimental sensitivity, /i —> ey, p —s- eee and muon conversion in a gold atom. The 


best current limits for the branching fractions are 54-56 


Br (/i —> ey) < 5.7 x 10 13 
Br (p —> eee ) < 1 x 10~ 12 


MEG 

SINDRUM 


(89) 

(90) 
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SINDRUM II. 


Br Au (fiN ->■ eN) < 7 x 10 


-13 


(91) 


The limit on muon conversion in gold is more stringent than current limits extracted 
using other nuclei, see e.g., 57,58] (titanium) or [59] (lead). 


In Section 2Y these three observables were expressed in terms of the Wilson coeffi¬ 
cients of the dimension-six SM effective theory Lagrangian, which in turn were deter¬ 
mined by integrating out the fifth dimension of the RS model. We now calculate the 
branching fractions of the observables for a given set of 5D input parameters. Before 
performing a scan through the parameter space of the model, it is useful to have some 
qualitative understanding of the effects of the various dimension-six Wilson coefficients. 
We will generally assume that the 5D Yukawa matrices are anarchic, without imposing 
any additional flavour symmetries. 


4.1 Estimates 

We first consider the effect of the dimension-six dipole operator, where we distinguish two 
different contributions: from the Higgs-exchange diagrams, which involve three Yukawa 
matrices, and from gauge-boson exchange, which involves only one. The gauge con¬ 
tribution leads to naturally suppressed flavour-violating couplings, whereas the Higgs 
contribution does not have a built-in flavour protection. For not too small Yukawa cou¬ 
plings the Higgs contribution is dominant. We mainly focus on /i —> e transitions, for 
which the dipole coefficients a \2 and « 2 i are relevant. 

To obtain an estimate of the Higgs-exchange contribution let us start with Wilson 
coefficient [see (69) and (80) 


a ij — 


Q„e r‘ 


■ /i“ > (l/r)[(2 + A^YY'Y + 


TO), 


(92) 


1927r 2 k 4 

in the custodially protected model with a bulk Higgs. An analogous expression holds 

For an exactly brane-localized Higgs, a^ is of 


(64) and (79 


in the minimal model 
similar size as above for the custoffially protected model, cf. (70), but suppressed in the 
minimal model due to the absence of the second line of (64). Now we recall that the 


relation of fermion zero-mode profiles, the 5D Yukawa matrix and SM Yukawa matrix 
(before rotation into the mass basis) is given by 


T 3 


= h 4h 1 / T ) i «9iffi/ T ) 


(o), 


(93) 


If the fermion mass hierarchy of the diagonalized SM Yukawa matrix is carried demo¬ 
cratically by left- and right-handed fermion modes, i.e. 


Vij 


y/m~i 


rru 


)/y/2 


we arrive at the estimate 


a H ~ 

*7 


192ir 2 t> 


[(2 + Ag K )Y? + B(; K y; 2 .], 


(94) 


(95) 
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where we assume that 


2 = [YY <Y]„ 

4 , - 


Yi 




2 = {y^Xk 

r u,* — v 

1 ij 


(96) 


are approximately independent of ij (’’anarchy”). For anarchic Yukawa matrices we also 
expect that the rotation matrices U and V follow the same hierarchy \U tJ \ m 
min ( y/rrii/mj , y/mj/rrii) and hence, barring accidental cancellations, it follows from ([ 9 ]) 
that a 12 = [U ] ]i k a k iVi 2 ~ a 12 . 

Further using that A^ K ~ 1.4 Bkk we obtain 14 


a 


Y?77. e 7?r^ 2 

1 n - I 

12 1927 t 2 v 


(97) 


which yields 


1 TeV 4 


Br (l-l ->• e 7) | Higgs dipole ~5-10 9 x — 


YJ 


(98) 


If the dipole also dominates ji —> 3e one can combine (13) and (15) to obtain the relation 

Br (/i —¥ 3e) 2 a e 


Br (n —y ey) 37 t 
which translates into an estimate of 


, rn M 11 

log —^ - — 
m P 8 


0.006, 


(99) 


Br (/i —> 3e) |Higgs dipole ~ 3 • 10 x 
For muon conversion one finds 


1 TeV 4 4 

J7 4 


Br A “ (/xJV -(■ eiV) |Hisg , dlpo]e ~ 1.5 x 10’ 11 x 


1 TeV 2 

J'i 


YJ 


( 100 ) 


( 101 ) 


We emphasize that these are crude estimates. Even in the anarchic case the random 
phases of the different elements can lead to cancellations or add coherently. However, 
they provide useful guidance to the results of the numerical scan discussed below. 

The Barr-Zee contribution is similar to the Higgs contribution, since the dominant 
contribution to the rj t j Wilson coefficient is also proportional to a product of three Yukawa 
factors. Comparing the prefactors in (37) we find that the Barr-Zee contribution to the 
dipole coefficient is smaller by a factor of about 170 than the contribution from the 5D 
Higgs loops. Thus we expect a /i —» ey branching fraction of about 


Br (/i ->■ ey) |BZ ~ 2 • 10 


7 3 1 TeV 4 

1 x Yf , 


4 


( 102 ) 


14 


Our estimates always yield a 21 ~ a 12 , hence we only give ol \2 explicitly. 
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Figure 13: Left: Gauge contribution to Br (// —> ey) (T = 1 TeV) for fixed Yukawa 
structure as a function of the absolute Yukawa coupling size. Right: The Wilson coef¬ 
ficient 7 ^ following the approximation (108) as a function of Yukawa coupling size for 
T = 1 TeV. The 0{rn il jrn T ') term is not included. 


if only the BZ contribution existed. The BZ contribution to the other processes is also 
smaller by a factor of about 170 2 . 

Due to the Yf dependence the Higgs-exchange induced dipole operator is less impor¬ 
tant for small Yukawa coupling. In this case, and also for the special case of the brane- 
localized Higgs in the minimal RS model, the dipole operator generated by gauge-boson 
exchange becomes crucial. We do not have an analytical expression for the gauge-boson 
contribution, but we know that there would be no flavour violation from it, if the func¬ 


tion Aij in (63) was independent of ij. The 5D mass parameters must decrease with 


the absolute values of the Yukawa couplings in order to guarantee the correct values for 
the SM masses fermion masses. Aij varies more strongly for smaller absolute values of 
the 5D mass parameters, see Figure [8j and therefore the flavour-changing gauge-boson 
contribution should increase with decreasing Yukawa coupling. To verify this we fix 
the Yukawa matrix structure, that is the ratios of all matrix elements, and scale the 
maximal entry Y max from 2 to 0.3. For simplicity we assumed symmetric 5D mass pa¬ 
rameters CL i = —Ce v The resulting /i —)■ ey branching fraction from oT alone in the 
minimal model is shown in Figure 13 (left). The precise value of Br (/r —>■ ey) obviously 


depends on the arbitrarily chosen Yukawa matrix structure, but the variation with the 
size F* of the Yukawa couplings is not very large compared to the fourth-power law of 


the Higgs-exchange contribution. For the Yukawa matrix used in Figure 13 we find a 
fx y ey branching fraction of a few x 1CT 12 . 

This agrees with the estimate based on the functional form of the gauge-boson induced 
The numerical value of the Wilson coefficient is 


dipole coefficient a?- 


15 


17 


v 


-6 (19) • 1CT 4 Uij 


(103) 


The value without (in) parenthesis is valid for the minimal (custodial) model and is 
independent of the details of the Higgs localization, y^ is the 4D Yukawa matrix in the 
flavour eigenbasis. The matrix relation a 9 cx y is only violated by corrections of about 
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(2-3)% as discussed in Section 3.3 This violation is the source of charged LFV as it 
introduces small off-diagonal elements in the dipole coefficients a tJ in the mass eigenbasis 


after EWSB. Using (94) and applying a factor 2/100 for the 2% of misalignment between 

i 9 

V?’ 


Dij and af-, we estimate 


ola ,12 ~ 2.6 (8.1) • 10 -8 x (104) 

for the coefficient relevant to fi —> e transitions. Again, we regard this as a rough 
estimate, since there may be cancellations when the rotation into the mass basis is 
performed. We then find: 


Br (/i ->• e 7 )| gauge dipole 
Br (/j -> eee) |gauge dipole 
Br Au (fJjN —>■ eA^) |gauge dipole 


r\_/ 






0.5(5) • 10“ n x 
0.3(3) • 10 “ 13 x 
0.2 ( 2 . 2 ) • 10" 13 


ITeV 4 

ITeV 4 

j'i 

ITeV 4 


x 


p4 


(105) 

(106) 
(107) 


Note that this contribution is independent of the typical size of anarchic Yukawa coupling 
up to the 0(1) variation shown in Figure 13 It is typically smaller than the Higgs 
contribution, but provides the “gauge-boson floor” to the dipole coefficient, since it is 
less sensitive to 5D model parameters than the Higgs contribution and always present. 
In the custodially protected model the rate is a factor of 10 larger than in the minimal 
model. 

The previous estimates were based on the assumption that the dipole operator dom¬ 
inates the LFV amplitudes. This is not always the case, especially for the yU —>■ 3e and 
muon conversion process. Next, we therefore consider the impact of the four-fermion and 
fermion-Higgs operators, which are generated at tree-level. In both cases the dimension- 
six Wilson coefficients are independent of the 5D Yukawa matrices. However, a depen¬ 
dence on the Yukawa matrices enters through the rotation to the mass basis after EWSB. 


For illustration we consider the operator (E^Ej) with Wilson coefficient 


c„ - = 


c}5ij and restrict ourselves to the minimal model 


For all three muon flavour-violating 
processes the relevant matrix elements are I4 2 - Flavour violation arises, because c} 
depends on the bulk mass parameter c Ei , hence c 4 - while diagonal is not proportional to 


the unit matrix in flavour space. We can estimate U^c^.V /2 by making use of hierarchical 
fermion zero-mode functions. Assuming f El (l/T) f E2 (l/T) -C f E3 (l/T) and symmet¬ 
ric mass parameters we can employ the rough estimate \V tJ \ ~ min (\Jrrii/rrij , yjmj/mi) 
with rrii being the SM lepton masses to obtain 


712 = Vlf) k V k2 ~ J— (cl -c\ + 0(m^/m T )). (108) 

Y 

We can use this formula to study the dependence of 7 ?- on the size of the 5D Yukawa 
couplings. Since the product of Yukawa matrix and 5D fermion zero-mode profiles must 
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reproduce the SM mass matrix to leading order in v/T, the 5D profiles and therefore 
the 5D mass parameters are correlated with the Yukawa matrix. The simplest estimate 
(assuming symmetric mass parameters) yields the correlation 1 /\fYl r^j where 

here Y* is the generic size of the anarchic Yukawa matrix element. In the following we 


do not distinguish this Y* from the one defined in (96) for order-of-magnitude estimates. 


Since the cf Wilson coefficients arise from a coordinate integral over a single fermion- 
gauge boson vertex they will roughly scale as [/^(l/T)] 2 , that is 1/Y*. This behaviour 
was already observed and explained in [13]. The right panel of Figure 13 shows 7 4 2 as 


-0.94 


a function of Y r (keeping the lepton masses fixed). The curve can be fitted by Y 
confirming the above scaling. This scaling is quite general, although if the mass hierarchy 


is mainly driven by the right-handed modes, the mass factors in the estimate (108) must 


change to account for the change in the relation (94). 


Similar estimates can be obtained for the four-fermion operator coefficients. Here 
we have terms with different dependencies on fiavour. The Wilson coefficient bf E of 


{L^^Li) (EffV'Ej) has three contributions denoted by b 0 , b\ and b 2 , see (43). b 0 does not 
depend on the 5D masses and hence does not contribute to flavour-changing processes. 
61 depends on a single bulk mass parameter and has the same scaling oc 1/Y* as the c a 
Wilson coefficients. The b 2 function depends on two bulk mass parameters and scales 
roughly as 1/Y 2 - However, as discussed below (52), for light leptons this term is sup¬ 


pressed and not relevant. This would be different for processes involving fermions with 
zero-modes that are IR brane localized such as the (right-handed) top quark, in which 
case the ratio of exponentials in b 2 no longer compensates the logarithmic enhancement 


factor and b 2 becomes the dominant term in (43). 


From (37), (22T) and (27 T) we see that the four-fermion coefficients usually appear 
in combination with the coefficients of the Higgs-lepton operators. For a typical RS 
model parameter point, which reproduces the lepton masses, the Higgs-lepton operator 
coefficients are larger by a factor loge relative to the four-fermion operator coefficients. 


This allows us to use (108) to estimate the effect of the tree-level operators on the 


generically tree-dominated LFV observables. We find 


„Au , 


1 — 12 


1-9 


ITeV 4 1 

(109) 

T 4 Y * 2 

ITeV 4 1 

( 110 ) 

T 4 Y * 2 ’ 


where we used the parameters given in Table [TJ We stress again that these numbers are 
rough estimates, which depend strongly on the precise structure of the flavour rotation 
matrices V and U. Thus we have three separate contributions to /i —> eee and muon 
conversion with different dependence on the size of the 5D anarchic Yukawa coupling 

«, U°. X' 2 )- 
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0.105658 GeV 

60 

s 2 

0.231 

60 

M z 

91.187 GeV 

60 

m t 

173 GeV 

60 

V[Au] 

0.189 

35 

S n [Au] 

0.0918 

35 

V" [Au] 

0.146 

35 

V[A\] 

0.0362 

35 

S n [Al] 

0.0167 

35 

V n [Al] 

0.0173 

35 

Jv P 

2 


fu 

Jv n 

1 


-CU 

Jp 

0.018 


£U 

J n 

0.016 


fs 

J n 

0.043 



M h 

M w 

r, 

S p [ Au] 
V p [Au] 

pAu 
1 capture 

S p { Al] 
V P [A1] 

pAl 

1 capture 

f d 

Jv p 

fd 

Jv n 

fd 
Jp 
fd 
J n 
fs 
J n 


5.10998 • 10“ 4 GeV 

60 


125.7 GeV 

60 


80.385 GeV 

60 


2.99598 • 10~ 19 GeV 

60 


0.0614 

35 


0.0974 

35 


8.71 • 10~ 18 GeV * 

6 

1 

6 

2 


0.0155 

0.0187 

4.64 • 10" 19 GeV 
1 
2 

0.034 

0.038 

0.043 


35 

35 

62 


Tabic 1: Input parameters for the numerical analysis. For the couplings of scalar quark 
currents to the nucleons we use the results of 36 and fix the value of the nucleon-pion 
cr-term to 50 MeV. The * indicates that we use the average of the values given in the 
references. 


4.2 Numerical analysis 

In the previous subsection we attempted to give an idea about the size and the relative 
importance of the various contributions to our three main LFV observables. While such 
estimates are useful to understand the rough dependence of our results on the input 
parameters, especially the Yukawa coupling size, they cannot replace a study of the full 
parameter dependence. To this end we next perform a numerical scan over the “generic” 
parameter space. We analyze four RS models: the minimal RS model (as defined in 115]) 
as well as a custodially protected model (as defined in |17|), each with either an exactly 
brane-localized or a bulk Higgs including its KK excitations in the (3 —» oo limit. We 
refer to these models as A4-bulk (minimal, bulk), A4-brane (minimal, exactly brane- 
localized), C-bulk (custodial, bulk) and C-brane (custodial, exactly brane-localized). 

The 5D input parameters needed for the numerical evaluation of the dimension-six 
Wilson coefficients are the 5D Yukawa matrices Y (and Y u in the custodially protected 
model), the 5D bulk mass parameters c^ = M^/k of the leptons and the KK scale T. 
In case of the exactly brane-localized Higgs the wrong-chirality Yukawa couplings can in 
principle differ from the “standard” correct-chirality Yukawa couplings, but for simplicity 
we assume them to be equal. The quark flavour parameters would affect our analysis 
only through suppressed terms, which are omitted (see Section [3]). 
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Since we do not want to give up the idea of “natural” Yukawa matrices, we further 
assume that the moduli of the complex Yukawa matrix entries are 0(1) and anarchic. To 
illustrate how the size of the 5D Yukawa matrix entries affects the different observables, 
we adopt two scan strategies. In both the modulus of the matrix elements is larger than 
0.1, but in the first (second) scan the maximal modulus Ym ax is bounded by 0.5 (3 for the 
second). The phases are arbitrary. Further, we require that the measured values of the 
charged lepton masses are reproduced by the chosen sets of 5D parameters. Contrary 
to the minimal RS model the custodially protected RS model includes a right-handed 
neutrino and a Dirac mass term for the neutrinos. Here we only require the neutrino 
masses to be below 0.1 eV. However, we do not demand that the PMNS matrix is 
reproduced by the 5D parameter sets, since the precise values of the neutrino masses do 
not affect charged LFV, and since, as mentioned in the introduction, the explanation of 
neutrinos masses and mixings is considered to lie outside the present model frameworks, 
as must obviously be the case for the minimal models. In practice, we randomly generate 
5D Yukawa matrices within the above mentioned constraints, and then fix the 5D mass 
parameters c^ t such that the correct lepton masses are obtained. For fixed value of the 
KK scale T and given scan strategy, we generate about 2-10 5 Yukawa matrices. For each 
of these we calculate the Wilson coefficients and then the branching fractions of fi —>• ey, 
H —> e conversion, /i —> 3e, r —* /jy, and r —* 3/i. The required low-energy parameters 
are shown in Table [lj We also added the material constants of aluminium, which serves 
as the target for the next generation of muon conversion experiments. 


4.2.1 Minimal model 


The results of our numerical scan through the constrained parameter space are best 
illustrated in two-dimensional scatter plots, which visualize the typical range of values 
for the branching fractions and correlations between the observables. It is important 
to keep in mind that the point densities in these scatter plots should not be used as a 
measure for the likelihood of the corresponding value in a given model. 


Figure 14 shows the values and correlation of the [i —ey and ^ —> 3e branching 
fractions in the minimal RS model for two different values of T, T = 4TeV (top) and 
T = 8TeV (bottom). T = 4TeV is also roughly the lower bound on the KK scale 
from electroweak precision observables 49 . The left panels correspond to the A4-bulk 


scenario, the right panels to the ARbrane case. Each plot shows the results for ax = 3 
in blue (dark grey) and for Y max = 1/2 in orange (light grey). The current and expected 
future experimental upper bounds are shown by solid and dashed lines, respectively. 

All four plots feature a sharp lower bound for /i —)■ 3e given the fi —> ey rate, which is 


precisely given by the relation (99). jd —> 3e branching fraction values in the vicinity of 


this bound are dominated by the contributions from the dipole operator. For very large 
dipole coefficients or equivalently very large /Me y branching fraction, the tri-lepton 
decay is always dominated by the dipole, and the two observables are strongly correlated. 
This generates the prominent thin line directed to the upper-right in the ARbulk model 
with large Yukawa couplings. 


40 







1(T 6 
1(T 8 
“ KT 10 
. ^- 12 
“ 1(T 14 
1(T 16 


I<4J 

i.10 


™1 

r M-brane, T = 8T< 

r 

r 

V 

1 ™1 

~l 

' "1 

Cl 

% 

1 

■■■■■ L... i 


^ ~\ 

r ' 

.j 


J__ 


10“ 


10 - 


10 12 10 - 
Br(/(->ey) 


10- 


10 - 


Figure 14: 2D scatter plots of the branching ratios // —> ey and p —> 3e for fixed 
T (4TeV upper row, 8TeV lower row) with Y max = 0.5 in orange (light grey) and 
Ymax = 3 in blue (dark grey), respectively. The left panels show the results for the A4- 
bulkmodels, the right panels for the AYbrane ones. The current experimental bounds on 
the branching fractions are given by solid lines. The region above/to the left is excluded. 
The sensitivity of future experiments is shown by the dashed lines. 


In the bulk Higgs case p —* ey is, as expected, quite sensitive to the upper bound F mav . 
This is a consequence of the YY'Y terms in the dipole coefficient. They are naturally 
flavour-violating and scale with F^ ax . Consequently, the scan with larger Yukawa entries 
includes points with substantially larger p —y ey branching fraction than the small 
Yukawa coupling scan. However, the dipole coefficient has two components. While 
the Higgs and the small Barr-Zee contributions scale as Y^ ax and vanish when the 5D 
Yukawa couplings go to zero, the gauge boson exchange contribution is not very sensitive 
to the Yukawa coupling size. In fact, for a generic anarchic Yukawa it grows mildly with 
decreasing Yukawa size, see Figure 13 Thus there has to be a smooth transition from 
the ’’Higgs-dominated” to the ’’gauge-dominated” regime when the Yukawa coupling 
decreases. 

To illustrate this point we included three curves in the plots for the bulk Higgs 
case defined as follows. We chose three (random) Yukawa matrices with Y max = 3 
and scale the matrices down to Y max = 0.25, keeping the relative size of the matrix 
entries fixed. The curves show the resulting trajectories. For large Yukawa couplings the 
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Figure 15: Correlation of the branching ratios for /r —>■ e y and p —> e (upper row) and 
for p —> 3e and p —> e (lower row), both for T = 4 TeV. The /e/t panels show the results 
for the Af-bulk, the right panels for the Af-brane model. 


curves all run close to the dipole dominance bound. With decreasing Yukawa couplings 
Br(/i — > ey) and Br(/i —>■ 3e) first also decrease following the change in the dipole 
coefficient. Then the growing effects of the tree-level operators begin to dominate p —> 3e 
and the corresponding branching fraction begins to increase, while Br(/i — > ey) continues 
to decrease. For even smaller Yukawa coupling the gauge-boson exchange contribution to 
Br(/i —* e'y) exceeds the rapidly decreasing Higgs contribution and Br(// —> ey) reaches 
a hard lower limit. 


The exactly brane-localized Higgs case displayed in the right panels of Figure 14 
behaves the same in this respect. However, since the leading Higgs contribution is 
suppressed for the exactly brane-localized Higgs, the range of values for Br(/i —> ey ) is 
almost independent of Y max . A change in Y m ax predominantly affects Br(p — > 3e), which 
increases for smaller due to the larger coefficients of four-fermion and fermion-Higgs 
operators. This also explains the drop shape of the scatter plot for Y max = 1/2. Points 
with large p —> 3e branching fraction arise from large tree-level Wilson coefficients either 
due to the structure of the Yukawa matrix or due to accidentally small couplings. In 
both cases the process p —> ey also receives sizeable contributions from the tree operators 
leading to the roughly linear correlation in the upper-right corner of the scatter points. 

Values and correlations of p —y 3e and p —> ey with muon conversion in gold are 
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shown in Figure 15 (colour coding as in the previous figure). The top row shows /i —y 
e 7 against fi —> e (Ad-bulk left, Ad-brane right). The two observables are essentially 
uncorrelated in the bulk Higgs scenario. This agrees with our previous observation that 
muon conversion is mostly insensitive to the dipole coefficient ad which governs the 
ji —> ey branching fraction. Only in rare cases is the dipole operator large enough 
to dominate also muon conversion leading to the noticeable spike towards the right 
in the upper-left plot. In the exactly brane-localized Higgs case (upper-right panel), 
correlations are absent only for Y* = 3. As mentioned before, for F* = 1/2 [i —> ey 
receives non-negligible contributions from tree-level operators, which manifests itself in 
a weak correlation. 

For the same reasons fi —» 3e and fi —> e (bottom row) are strongly correlated 
for small Yukawa couplings in Ad-brane model, but only feature a lower bound on the 
branching fraction of [i —> 3e for a given Br(/i —y e ) in the other scenarios. As noted in 
the previous subsection, the branching fraction of /r —> e decreases with increasing values 
for Y*. This effect can best be seen in the upper left panel of Figure 15 The slopes of 
the three sample trajectories also verify this effect. 


4.2.2 Custodially protected model 


Figure [16] shows the combined results for the custodially protected model. The left 
panels correspond to the bulk Higgs model C-bulk and the right panels to the exactly 
brane-localized model C-brane. The colour coding is the same as above. Here the KK 
scale T was fixed to 8TeV, since for T around 4TeV it is already non-trivial to find 
points, which are not in conflict with the muon conversion bound. 

The broad picture for model C-bulk is almost the same as for Ad-bulk. The shape 
of the distributions does not change. Quantitatively, the custodially protected model 
generates significantly larger branching fractions. In particular, the /r —>■ ey branching 
fraction, which is most sensitive to the magnitude of the dipole operator coefficient, 
is typically enhanced by a factor of about five. This was to be expected as the main 
difference in the custodially protected model is a larger gauge- and Higgs-contribution 
to the dipole coefficient ad. Again the dipole operator creates a correlation of (i — > ey 
and /x ^ 3e especially for the larger value of fi iaY . 

For the C-brane model the distributions are very different from the result in the 
Ad-brane model. This is a consequence of the new additional terms (70) in a fi This 


additional contribution to the dipole coefficient is only slightly larger than the corre¬ 
sponding contribution for a bulk Higgs. The phenomenology for bulk and brane Higgs 
case is therefore quite similar in the custodial RS model. The bounds imposed by the 
non-observation of LFV are comparable, although more restrictive for the exactly brane- 
localized Higgs. 

The fact that the sign of a fi depends on the Higgs localization does not lead to a 
noticeable effect. If the dipole operator is dominated by the Higgs contribution, a sign 
flip of the coefficients a^- only affects terms in (13), (15) and © that come from an 
interference of the dipole with a four-fermion operator. In general, these terms do not 
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Figure 16: Correlation of the branching ratios for /i e 7 and fi —> 3e (top row), 
H e 7 and fi —> e (middle row), and n —> 3e and /a —>■ e (bottom row) for T = 8 TeV in 
the custodially protected model. Model points with Y max = 1/2 are indicated by orange 
(light grey) points. For Y max = 3 we use blue (dark grey) points. The left panels show 
the results for the C-bulk, the right panels for the C-brane model. 
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provide the dominant contribution to the branching fractions. The situation would be 
different if the RS contribution could interfere with a sizeable SM contribution to LFV 
observables. An observed enhanced or reduced rate could then be used to discriminate 
the brane from the bulk Higgs model p*] 


4.2.3 EDM constraint 


The randomly sampled Yukawa matrices also generate electric dipole moments (EDMs) 
of the leptons through the non-hermitian part of a-, see (|8]) . In case of the electron the 
relation is given by 

d e = F 3 (0) = m e i(A R -A L ), (111) 

where the form factor now refers to the flavour diagonal electron-photon vertex, and 
A l r are defined as in (37), (38) with subscripts 12 replaced by 11, and m M by m e . 


We checked that the present experimental limit on the electron electric dipole mo¬ 
ment 1631, 


\d e \ < 8.7 -10 29 e cm (at 90% CL), 


( 112 ) 


does not affect our conclusions. That is, while up to 90% of the randomly scanned 
model points in the sample with Y nia Y = 3 (and up to 25% in the Y max = 0.5 sample) 
fail the EDM constraint for the bulk Higgs models, the EDM and the LFV observables 
are uncorrelated, so that the ranges covered by the scatter plots look almost exactly 
the same, when the points failing the EDM constraint are excluded. Only the tips of 
some of the spikes pointing to the upper-right in the scatter plots are cut-off in the 
Y max = 3 sample, but these points are also excluded by LFV constraints (see solid lines 
in the previous figures). The predicted values of \d e \ and the /i —> ey branching fraction 


are shown in Figure 17 for the different RS models and values of T considered in this 
analysis. The featureless shape of the area filled by the sampled points illustrates the lack 
of specific correlations. The scaling of both observables with Y* is visible (orange/light 
grey vs. blue/dark grey sample) whenever there is an unsuppressed contribution from 
Higgs exchange. 

We can turn this question around and ask what are the expectations for value of 
the electron EDM after the LFV constraints on the RS models are taken into account. 
To this end we create histograms of \d e \ for our samples before and after applying the 
present LFV constraints. The Ad-brane models are special, since the Higgs contribution 
is strongly suppressed. In these models the EDM distribution peaks at a few times 
10 _3 °e cm for T = 4TeV, well below the present limit (112) and is hardly altered by 
the LFV constraints. The tail of the distribution above about a few times 10 -29 e cm 
is, however, cut away independent of the value of the Yukawa coupling. This is different 
for the other models where the distribution and impact of constraints is sensitive to the 
maximally allowed Yukawa coupling Vnax- This is illustrated on the example of the Ad- 
bulk model in Figure [T8j where the left two panels refer to T = 4TeV, Ym ^ = 0.5 and the 


15 


This is precisely what is observed in Higgs production, see e.g. 131 . 
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Figure 17: Correlation of the electron EDM with the /i —)■ e 7 branching ratio for all 
six scenarios. Data sets with Y* = 1/2 are indicated by orange (light grey) points. For 
17 = 3 we use blue (dark grey) points. The left panels show the results for the bulk 
Higgs, the right panels for the exactly brane-localized Higgs. 
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Figure 18: Histograms of the electron EDM distribution of the A4-bulk sample before 
(orange/light grey) and after (darkgrey/blue) the LFV constraints are applied for T = 
4TeV, Ymax = 0.5 (left) and T = 8TeV, Y max = 3 (right). The first row shows the 
absolute number of representatives of the sample, the second the normalized distribution. 
The vertical line indicates the current bound (112) on the EDM. 


right ones to T = 8TeV, T max = 3. The respective upper plot shows the distribution of 
d e before (light grey/orange) and after (dark grey/blue) the LFV constraints are applied. 
The second distribution is hardly visible, as most of the models in these two samples are 
excluded by the LFV constraints. In the second row we therefore show the normalized 
distribution for the two cases. In the first case, T = 4TeV, V max = 0.5 the tail of larger 
values of d e is cut away, while the peak of the distribution remains at the same value 
as without the LFV constraints. In the second case, T — 8 TeV, V max = 3, most of 
the distribution without LFV constraints lies above the EDM bound (112) indicated by 
the vertical line. When the LFV constraints are taken into account the expected value 
of the electron EDM reduces by more than an order of magnitude. We note that with 
LFV constraints applied the current EDM bound is always close to the upper end of the 
distribution. In other words, quite generally, the lepton-flavour violating processes put 
similar generic bounds on the electron EDM as the current bound from the direct EDM 
measurement. Of course, these conclusions hold in the parameter space of RS models 
with anarchic Yukawa couplings in a statistical sense, and not for any particular model. 
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Figure 19: Correlation of the branching ratios of r —> /17 and r —> 3/i for T = 4TeV 
in the minimal model. The left panel shows the bulk Higgs, the right the exactly brane- 
localized Higgs case. 


4.3 A note on LFV r decays 


Tau decays offer another opportunity to study LFV. However, the short lifetime of the r 
and its high mass make it unsuited for studies in low-energy facilities. The best bounds 
on processes like r —> e^y or r —>• 3/i come from Babar 
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Belle 165] as well as LHCb 1 66 

The RS model naturally generates higher rates for r — > fi, e transitions than for 
fi —> e transitions, since there is a close relation of lepton masses with the corresponding 
zero-mode profiles, which also control the size of LFV. However, the fantastic sensitivity 
of past and future experiments searching for muon flavour violation still makes searches 
in the muon sector the most promising avenue, unless an additional flavour structure 
suppresses muon flavour violation. 

Nevertheless, it is instructive to provide the expectations for tau LFV in the RS 
model. In Figure 19 we show the values and correlation of the r —> /zy and r —> 3/z 
branching fractions. The colour coding is the same as in the previous subsection. The 
solid lines correspond to the current best upper bounds on the branching fractions. 
Compared to the bounds in the muon sector the current limits from tau decays are not 
restrictive even for T = 4 TeV. An improvement of more than five orders of magnitude 
would be required for constraints as severe as those from muon decays. 

Qualitatively, the r —» /17 vs. r —» 3/i plot is similar to the corresponding “muonic” 
plots (first row in Figure |l4|). The main difference is the large effect of the four-fermion 
and fcrmion-Higgs operators. For the exactly brane-localized Higgs this generates the 
strong correlation for small Yukawa couplings Y max = 0.5. In the bulk Higgs case this 
effect prevents scatter points close to the dipole-dominance line. 


4.4 Discussion 

In the following we summarize and emphasize the main conclusions from the phe¬ 
nomenological study. Based on the topologies of the decays we (naively) expected 
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that fi —» e conversion in nuclei and fi —> 3e are driven by the current-current oper¬ 
ators (#>), whose Wilson coefficients are determined by the tree-level four- 

fermion and Higgs-fermion coefficients, c“- and bij, in the dimension-six SM effective 
Lagrangian ([3]). Whereas the rate for radiative muon decay should primarily be deter¬ 
mined by the loop-induced dipole operator. This situation would be ideal, as the two 
groups of processes would then provide complementary information on the underlying RS 
model through the different parameter dependencies of the relevant Wilson coefficients. 

However, and this is our main message #1, we find that this is not generally the case. 
In particular, the process fi —>■ 3e can receive sizeable contributions from dipole operators. 
This can even be the dominant effect as indicated by the prominent dipole dominance 
line in the plots. In the minimal RS model with an exactly brane-localized Higgs and 
small Yukawa couplings the situation can also be reversed. Then fi —> ey receives relevant 
contributions from tree-level coefficients and is no longer governed by the dipole operator 
alone. This makes it clear that it is necessary to consider all dimension-six operators for 
every observable to arrive at a reliable picture of LFV in the RS model. 

From our numerical analysis we also see that it is important to include both, the Higgs 
and the gauge-boson exchange contribution to the penguin diagrams that generate the 
dipole Wilson coefficient a because they exhibit a fundamentally different dependence 
on the parameters—most notably the Yukawa coupling. Only the full dipole coefficient 
including all contributions gives an accurate description of the fi —> ey decay over a 
large range of Yukawa sizes. Without the Higgs contribution the dominant effect for 
large Yukawa couplings is completely missed. Whereas the presence of the gauge-boson 
exchange contribution prevents the dipole from becoming irrelevant for small Yukawa 
couplings. This manifests itself in a lower bound on Br(/r —)■ ey), which depends only 
on the structure of the Yukawa matrix but not on the overall size of the couplings. This 
is our main message #2: in RS models Higgs and gauge-boson induced flavour violation 
can both be important and only a full calculation of the dipole operator coefficient gives 
a reliable picture of LFV over all of the parameter space. 

The scans suggest that anarchic RS models with minimal particle content typically 
need KK scales T larger than 4 TeV to be compatible with the current data on charged 
LFV. For models with a bulk Higgs the combination of fi —> ey and muon conversion 
makes the bound on T almost insensitive to the size of the Yukawa matrix. For an 
exactly brane-localized Higgs the situation is more complicated. T > 4 TeV is still valid 
for small Yukawa couplings. However, in this model the limit on T comes mainly from 
muon conversion data. As discussed above, muon conversion is dominated by tree-level 
operators whose Wilson coefficients decrease for increasing Y*. Thus, the bound on T 
gets weaker the larger the Yukawa couplings. For Y max = 3 we can still fold numerous 
data point that satisfy all constraints for T as low as 2 TeV. In the custodially protected 
model on the other hand, the larger particle content generally leads to larger Wilson 
coefficients. Consequently, the lower bound on admissible KK scales T is higher. It 
is noteworthy that in the custodially protected model the bound on T is essentially 
independent of the Higgs localization for all three observables. 

The limits on the KK scale T will significantly improve in future experiments. Espe- 
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daily the next generation muon conversion searches will provide strong constraints on RS 
models. Our analysis was performed for a gold target nucleus 
and COMET 
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Mu2E 68 


-DeeMe 

use silicon and aluminium target nuclei, respectively. Aluminium and 
silicon have an approximately 20 times smaller muon capture rate compared to gold, 
which enhances the branching fraction, but the wave-function overlap integrals (T>,V,S) 
relevant to muon conversion are also smaller, see Table [lj Thus an expected lower bound 
on the branching fraction of 7 x 10 -17 for aluminium is roughly equal to a bound around 
10 -16 in gold 1351 (indicated by the dashed lines in the plots). This combined with the 


expected improvement on Br(/i —» ey) during the next run of the MEG experiment 70 


could exclude the parameter space of anarchic RS models (custodially protected or min¬ 
imal) up to a lowest KK resonance mass of 20 TeV, which corresponds to T > 8 TeV. 
The current /i —>■ 3e constraint is less constraining than the one from muon conversion. 
However, the proposed Mu3e experiment 71 aims for sensitivity to a branching fraction 


of about 10~ 16 . At this level /i —>■ 3e alone will be able to exclude anarchic models with 
T < 5 TeV. Naturally, should LFV be observed in any of the new experiments, the 
model-dependent correlations among different processes can be used to further constrain 
the RS parameter space. Hence, our main message ^3: LFV violation provides very 
strong constraints on RS models, and future experiments will further strengthen them. 

ft is interesting to compare the charged LFV constraints on the RS model and its 
KK scale to those derived from other processes. The non-observation of direct KK gluon 
production forces T to be larger than only about 1 TeV. This cannot compete with the 
bounds from electroweak precision observables, notably the S and T parameter, which 
are essentially model-independent. They only depend on the particle content of the 
model and to a lesser degree on the 5D Higgs profile 72 . For the two models discussed 


in this work the electroweak precision observable bounds are T > 2.3 TeV (custodially 
protected) and T > 4 TeV (minimal) [49j. For not too large Yukawa couplings the RS 
contribution to (g — 2) M is model-independent in the same sense 15 17 . However, in 


this case the SM contribution is non-zero and given the present situation of experiment 
and SM theory, the bounds on the KK scale are not competitive. 

If one allows for a somewhat stronger dependence of the bounds on the model pa¬ 
rameters, Higgs production (and subsequent decay) is also an interesting observable, ft 
depends more strongly on the (mainly quark) Yukawa matrices Y u ,Y d than the oblique 
parameters, but is still far less sensitive to its detailed structure than processes like 
/i —> 3e, because Higgs production depends to leading order only on the traces of Y^Y U , 
Y r j Y d . The trace of a product of anarchic matrices follows a narrower distribution than 
an individual matrix element. One Ends that T has to be larger than 2 (4) TeV @ 95% CL 
for Y ~ 3 in the minimal (custodially protected) model with a narrow bulk Higgs |3Tj. 
For smaller Yukawa couplings the bound becomes weaker as the effect on the production 
cross section decreases with Y*. For the exactly brane-localized Higgs the constraints are 
stronger, and one finds the same bounds on T as above already for T* m 1. As we have 
seen the situation is different for LFV observables. In the minimal model an exactly 
brane-localized Higgs leads to weaker bounds than a bulk Higgs, and in the custodially 
protected model the bound from the bulk Higgs (with KK modes) is comparable to 
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the one in the exactly brane-localized scenario. Thus for large quark Yukawa couplings 
and the exactly brane-localized Higgs, Higgs production provides at least equally strong 
bounds on the KK scale than the non-observation of charged LFV. The LHC will be 
able to improve on this further in the future. In all other cases the next generation LFV 
experiments will be able to set the most stringent limits on the KK scale. Of course, 
this comparison assumes that the magnitude of anarchic Yukawa couplings is roughly 
the same in the quark and lepton sectors. 

The quark-sector Yukawa couplings also enter the RS modification of meson oscil¬ 
lations, where in particular ex, which measures CP violation in kaon mixing, is very 
sensitive to coloured KK states [8.,:9j. If we follow (8| and estimate this effect from the 
dominant left-right four-quark operator ( Qi'y^T A Qi){D j^T a D j ) along the lines of Sec¬ 


tion |4.1[ we find that T in excess of 8 TeV is needed to avoid conflict with experimental 
data. We have seen in the previous section that for any given data point /i —> ey can 
differ by orders of magnitude from the simple estimate, since cancellations may or may 
not be present for a given set of 5D Yukawa and mass parameters. The same is true 
in the quark sector and much smaller modifications of ex than indicated by the naive 
estimate are possible, see |9). We confirmed this independently using the results of |73|. 
Thus, from an analysis similar to the one performed in the previous section, as well as 
from (9j, the limit of the TeV brane scale from ex is around T = 7 TeV — quite close to 
the one expected from the next generation LFV experiments. However, contrary to ex, 
whose ultimate reach is limited by the theoretical precision of the SM prediction, the 
potential of LFV observables, in particular of /1 —)■ e conversion, is not yet exhausted. In 
any case, given that the quark and lepton Yukawa couplings may exhibit different pat¬ 
terns, independent tests of the RS model in both the quark and lepton flavour sectors 
in a wide range of observables are not redundant. 


Our findings can be compared to the results of 113 , which provided the first detailed 


analysis of lepton flavour violation in the minimal RS model in the KK mode picture. 
The branching ratio of /i —> ey is determined from the Higgs-exchange contribution to 
the dipole Wilson coefficient alone, which is computed via one-loop diagrams involving 
the Higgs zero-mode and first fermion KK excitation. Muon conversion and /i —> 3e 
are computed from the tree-level Wilson coefficients, while the dipole contribution is 
neglected. Both exactly brane-localized and bulk Higgs scenarios are investigated. We 
can compare with our full results only for the bulk Higgs case as their exactly brane- 
localized Higgs result for the dipole operator is cut-off dependent. Despite these caveats 
the overall size of the bound on the KK scale T for a bulk Higgs is compatible with the 
one found above. The main difference is the dependence of the branching fractions on the 
model parameters. As 13 only includes the Higgs contribution J^j the dipole coefficient 
has a straightforward dependence on the Yukawa coupling size Y*. The identification 
of the dependence on the Yukawa coupling size as a distinguishing feature of tree-level 
and loop-induced observables, i.e. muon conversion and tri-lepton decay as opposed to 
t’ —» f'y, is however valid only for medium-size Yukawa couplings, since otherwise the 


16 


The absence of the KK modes does affect the qualitative characteristics of the Higgs contribution. 
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neglected gauge-boson contribution with its different dependence on Y* becomes relevant 
for small Y*, and for large Y* the tri-lepton decay is dominated by the dipole operator 
and therefore effectively also loop-induced. 

We can also compare our results with [14]. Here the Higgs-exchange contribution 
to the dipole operator is not considered. The dipole coefficient is computed in the 5D 
framework from a subset of gauge-boson exchange diagrams including a dimension-eight 
effect with three Yukawa matrices. Comparing orders of magnitudes their results for 
/j, —> ey and /i —)• 3e are similar to our exactly brane-localized Higgs case in the minimal 
model. In particular, the lower bound on the fi —> ey branching fraction for small 
Yukawa couplings is also present in their estimates. 


5 Summary 


In this paper we have undertaken a comprehensive study of charged lepton flavour vio¬ 


lation in Randall-Sundrum models. The approach follows the strategy developed in 15 


that is we assume that the KK scale is significantly larger than the electroweak scale and 
match the RS theory to the SM effective theory including dimension-six operators by 
integrating out the extra dimension in a fully 5D quantum-field-theoretical framework. 


Some of the Wilson coefficients of the dimension-six operators could be taken from 15,17 


and the remaining ones were computed here. We considered the RS model with minimal 
field content and an extended model with additional fields to protect low-energy precision 
measurements from custodial symmetry violating contributions. We further considered 
three implementations of the IR-brane localized Higgs field. The exactly localized case 
and the limiting cases of a bulk Higgs with and without KK Higgs excitations. The third 
scenario was motivated by the recent work 20 , which demonstrated the non-decoupling 


of Higgs KK modes in the KK mode picture. We confirm this finding within the 5D 
framework. 

Our calculation is the first complete one of dimension-six effects (more precisely 
those not suppressed by powers of small lepton masses), and considerably sharpens 
previous results from [l3 ) 14,[T6|. This concerns in particular the Wilson coefficient of 
the electromagnetic dipole operator L l a lw EjF lw . Not only does it depend on the way the 
Higgs is localized near the IR brane, which has phenomenological consequences. It also 
receives three contributions with different dependence on the magnitude of the anarchic 
5D Yukawa matrices, which can all be important in certain parameter regions. Amongst 
these the gauge-boson exchange contribution is computationally the most demanding. As 
already emphasized in the discussion of our results in the previous section the interplay 
of the three contributions leads to distinctive features in the scan of the parameter 
space and it is important to include them all. We also find that /i —> 3e can receive 
sizeable contributions from dipole operator, while in some cases /i —> eq receives relevant 
contributions from tree-level operators and is no longer governed by the dipole operator 
alone. 

Assuming generic anarchic Yukawa matrices we studied the typical range for the 
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branching fractions of fi —> ey, /i —> 3e, jiN —> eiV as well as r —)■ / 17 , r —>■ 3/i in both 

the minimal and the custodially protected RS model. The combination of /i —)■ ey and 
gN —)• eiV currently provides the most stringent constraints on the parameter space of 
the models. A typical lower limit on the KK scale T is around 2 TeV in the minimal model 
(up to 4 TeV in the bulk Higgs case with large Yukawa couplings), and around 4TeV in 
the custodially protected model, which corresponds to a mass of up to 10 TeV for the first 
KK excitations, far beyond the lower limit from the non-observation of direct production 
at the LHC. The next-generation LFV experiments will push the lower limit on the KK 
scale T further up. Given their projected sensitivity each of g —* ey, /iN —> eN and 
/i —>■ 3e will contribute to improving the current bound. When combining all searches, 
the non-observation of lepton-flavour violation will exclude anarchic RS models without 
additional flavour symmetries up to a KK scale of T ~ 8 TeV, which corresponds to KK 
gluon masses of about 20 TeV. 

We also correlated the electric dipole moment predicted by the RS models with lep¬ 
ton flavour violating observables in the presently considered scenario of random anarchic 
Yukawa matrices. We find that the non-observation of charged LFV in current experi¬ 
ments imposes a similar bound on the electron EDM as the one set by the direct EDM 
measurement. 
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A The bulk Higgs 


For the bulk Higgs we follow |47,48]. The 5D Higgs action reads 


So, — d 4 x / dz 


(kz) 


g MN {D M ^ {D N <S>) - 


3(\/d55( z — 1 /T))Vi/t — d(y/g55(z — l/k))Vi/k 


(113) 


V\/t — —7711 /t T A 


The brane potentials are 

Vl/k = 777i /fc 4> t 4>, 
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(114) 








with rrii/k = (2 + j3)k. We define /3 = \J A + /i 2 . Note that 
is, up to normalization, written as 


m 


47 


hi/T — 


A 

2P 




y TeV 


the IR brane potential 


(115) 


where the coupling constant A is dimensionless. For our purposes the form of (114) 


is more convenient. The choice for the UV potential parameter rrii/k leads to a Higgs 
vacuum expectation value (vev) that rises towards the IR brane for positive j3. 

A.l Zero mode and vacuum expectation value 

The zero-mode equations of motion are given by 


3 d z z~ 3 d z - <3> (0) (z) = -m 2 0 &°\z) 


dM 0) 


$( 0 ) 


z=l/k 


/k •> 


a $(°) 


$(o) 


z=l/T 


= m i/T 


T 

k 


(116) 

(117) 


The boundary conditions ensure that the boundary terms arising from integration by 
parts vanish, is the squared mass of the zero mode in the unbroken phase. As we will 
see below it is tachyonic and of the order of the physical Higgs mass, i.e. much smaller 
than the KK scale T. 

The general solution of the differential equation is 


$ (0) (z) = A f 0 z 2 (Jp(m 0 z) + Cy 0 (m o z)) . 
The UV-brane boundary condition can be used to determine 


C = 


Jp+i(m 0 /k) 


2- 2 ( 1+ ^7T /mo \ 2+2/3 


for ?n 0 k . 


yp+i(mo/k ) T(l + /3)T(2 + /3) V k ) 

Observing that Jp(x)/yp(x) oc x 2/3 for small arguments, we find that 

J/ 3 (rn 0 z) > Cy 0 (m o z) 

for mo T. We can use this approximation for the zero mode to obtain 
<h (0) (^) = Af 0 z 2 Jp(m 0 z) fa A/"o z 2+l3 + O(m 0 z ). 

The overall normalization is given by 




(118) 


(119) 


( 120 ) 


( 121 ) 


( 122 ) 
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Up to higher terms in uiq/T the zero-mode mass is determined by the equation 




<f> 


= T {2 + (3) - m 0 


z=l/T 


Jp+\{mo/T) 

Jp(mo/T) 


T 


my T - 


Expanding the Bessel function for small argument, we find 


m i /rj - T(2 + P) = 




2(1 + /3)T 


+ higher-order terms, 


which implies 


? 7 ln ~ 


2 (i+« ( 2 ++- 


(123) 


(124) 


(125) 


Note that mo must be small compared to T, otherwise the expansions above would not 
have been allowed. We return to this point below. 

The 5D profile of the vev is not needed in our computation, since it is done in 
the unbroken electroweak phase. The Higgs vev only enters at the 4D level in the 
effective Lagrangian—as a low-energy parameter determined from experiment. Still it is 
instructive to see how the vev profile arises. To this end we substitute <f> —* ~j^(v + h) 
into the Lagrangian (113) and expand all terms (see 31 for a more detailed derivation). 
We can use that 4D derivatives on v vanish. This leads to the equations 


T T 

{-d z + —nrii/T - 3-\v 2 )h\ z ^ 1/T = 0 

(9z TTl\/k)h\z—¥\/k 0 

T T 

(~d z + - —\v 2 )v\ z ^ 1/t = 0 


(126) 

(127) 

(128) 
(129) 


(dz ~ mi/k)v\ z ^i/k = 0 
along with the standard equation for Higgs bulk profiles. This gives the solution 

v(z) = J\f v z 2+13 , (130) 


which is strongly IR localized already for moderately large, positive values of (5. The IR 
boundary condition determines 

K 2 = (mi/T - (2 + P)k) ^T 4+2/3 • (131) 

Equivalently, by requiring that W boson acquires the correct mass 

AC = v'rvpiw M. (132) 

With this input we can compute the physical Higgs mass 

m H — m l + 6(1 + /3)^ttXv(1/T) 2 

K 


55 











(133) 


= 2(1 + 0) ((2 + P)k - m 1/T + 3At»(l/T) 2 ) . 

Using 


j! (l/Tf = m s _ (2+ffl 


k 


this result can be rewritten into 


k 


v(l/Tf 


m), = 4(1 + P)\ ' ’ T 2 « (125GeV) 2 


and 


m 


H 


(134) 


(135) 


(136) 


Thus we find \rri^\ -C T 2 , which was necessary to justify the expansion in the broken 
phase. We note that the requirement that mu ~ 125 GeV implies a fine-tuning between 


the parameters nri\/T and (2 + f3)k, see (125). We further note the relations 


Au(l/T) 2 = 


m 


H 


4(1 + (3)T 2 


k, 



mi/ T = (2 + /3)k + 


A.2 Higgs Propagator 

The 5D Higgs propagator is determined by the equations 


(137) 

(138) 


P 2 ~ ^ + z 3 d z z~ 3 d z 


A$(p, z, z') = i{kz') 3 5(z — z '), 


d z A®(p,z, z')\ z = l/k = mi/kA$(p, 1 /k,z'), 

T 

d z A t (p,z,z)\ z=1/T = mi/r-j- A*(p, 1/T, z ), 


(139) 

(140) 

(141) 


which can be solved in the standard way, see e.g. 15 . After Wick rotation to euclidean 
space, the full Higgs propagator is given by 

A ^{p, z, z’) = Q(z - z’) ik 3 z 2 z ' 2 

(P+i ® Kci(pz’) + Kp+i (£) Ip(pz')) (I „ +1 (p Kefpz) + A > +1 (f) Ip(pz)) 


X 


fft. (!) Kp +1 (p - Kp +1 (|) I „ +1 (p 


+ {z -f->- z'} , 


(142) 
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where K and / are modified Bessel functions. It is useful to not only have the full 
propagator, but also the zero-mode subtracted propagator. We only work to leading 
accuracy in v/T, that is we approximate 

m l/T j^ = (2 + f3)T . (143) 

The Higgs zero mode is then massless, and its profile is proportional to the vev profile 
derived previously. The zero mode can readily be removed from Euclidean propagator 
via 

Af s (p, z,z') = A*(p ,z,z<) - (144) 

(~P ) 

since removing the zero mode corresponds to removing the pole at p 2 = 0 from the full 
propagator. 


A.3 Yukawa matrix scaling 


For the bulk Higgs field the Yukawa coupling develops a dependence on the Higgs 5D 
mass /i or, equivalently, j3. To see how this dependence arises let us compare the situation 
with the delta-regularized narrow bulk Higgs (40). In the latter case, we find for the 4D 


SM lepton Yukawa matrix the standard expression 


Vij = fL-( 1 / T )9E’( 1 / T )l^ Y ij + hi £ her terms 


,( 0 ) 


T 3 


fc 4 


1 - 2c Li 


e 1 2cL i \l 1-e 


1 + 2ce 

m y. 

1+2 c E . 1 ij ■ 


(145) 


For the bulk Higgs the bulk action contains the interaction term 

where L, E, $ are 5D fields, and Y 13 is the dimensionful bulk Higgs Yukawa coupling. 
Inserting zero modes and integrating over z, we obtain (up to terms suppressed by powers 
of e) 


drxYij Li{x , z)Q(x, z)Ej{x , z) + h.c. 


(146) 


Vij 


1 — 2 CLi I 1 + 2cEj a / 2(1 + ft) k 1 / 2 
-; 


(147) 


Since the SM Yukawa coupling should remain finite for large /3, the bulk-Higgs Yukawa 
coupling scales as 


y/3 2 “ C Li + C Ej + P P~*oo VP 

1 x VWW) V2' 


(148) 
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Figure 20: Left panel: Diagram contributing to the matching onto L{Ej L k Ei. The inter¬ 
mediate Higgs propagator is zero-mode subtracted to remove long-distance contributions. 
Right panel: Wilson coefficient Cf{ E f E as a function of f3 for T = 1 TeV, cl 1 = — Ce x = 0.6 
and Y — 1. 


Comparing the expressions (145) and (147) we identify 


Y P = Y 2 ~ C E + Ce j + P 
y/k 2(1 + /?) 


(149) 


A.4 KK Higgs example: Four-fermion operators 

To gain some intuition for the properties of the Higgs KK modes we consider the example 
of the Feynman diagram in Figure [20j which might contribute to the matching of four- 
fermion operators of the form L^Ej L k Ei. The corresponding Wilson coefficient is given 

by 


C m LE _ dx [T dy ( 0 ) 


_— iY p Y^ . _ , _ 

T 2 ij kl J i (kxfji (kyf 


/ff(*)^(*)^T(p = ^x,y)f^{y)g^{y) . 


, ZMS 


(°V„,U 0 b 


(150) 


For vanishing four-momentum exchange the zero-mode subtracted Higgs propagator has 
the particularly simple form 


A™ s (0 ,x,y) = 


ik 3 x 2 Py 2 P 

w 


TW b ( T 2 (* 2 + y 2 ) - 2) - ^ 


+ x 2l3 0(y — x) + y 2l3 d(x — y) 


(151) 


where we dropped terms suppressed by powers of T/k. With this expression the integrals 
over x and y in (150) are straightforward, and Ck EEE can be determined analytically for 


all values of /3. C^ EE vanishes as 1//3 in the limit /3 —» oo, as illustrated in the right 


panel of Figure 20 
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This result can be understood by looking at the defining expression. The scaling with 
(3 is determined by three factors: the Yukawa matrix scaling, the scaling of the Higgs 
propagator, and the scaling of the integration variables x, y in the relevant integration 
regions. The two Yukawa couplings each contribute a factor of \f]3. The Higgs propa¬ 
gator is slightly more complicated. Let us examine the three terms square brackets in 
(151) separately. The first term (without step-functions) does not feature an immediate 
suppression for large j3, since the 1 / f3 in the prefactor of the square bracket is cancelled. 
The suppression arises only after integration over the bulk coordinates. To see this, 
we write x and y in the overall factor ( T 2 xyY as 1/T(1 — x 0 / (3) and 1/T(1 — yo//3), 
respectively, such that Xo and yo measure the distance of x, y from the IR brane in units 
of 1/((3T), the typical scale for Higgs KK excitations. We then find factors of the form 
(1 — xq/ fiY and (1 — yof/3) 13 , which behave as e~ x ° and e~ vo for large (3, respectively. 
Hence the first term counts as 0(1) only for x and y within 1 /((3T) of the IR brane. 
The 5D coordinate integrals then count as l/(/3T) each, and the contribution of the first 
term in (151) to the Wilson coefficient is of order (\//3) 2 x 1 x 1 //3 2 , which vanishes for 
large (3. The remaining two terms in (151) have different properties. There is a global 
factor of 1//3, but there is no requirement that x, y are close to the IR brane. Let us 
focus on the second term, which is non-zero only for y > x. It contains the factor ( x/yY, 
which ensures that the contribution to the Wilson coefficient is exponentially suppressed 
if x -C y (1 — 1/(3). Changing integration variables from x, y to y, x — y shows that the 
integral over x — y counts as 1/((3T), while the integral over y is effectively unconstrained. 
Thus the overall total scaling is (v(d) 2 x 1//3 x 1//3, which also vanishes for large f3. The 
same argument with xfti/ ensures that the third term in (151) does not contribute to 


B KK Higgs contributions 

The contribution of KK Higgs modes to a^ is proportional to the corresponding contri¬ 
bution for the Higgs zero-mode for each diagram topology. It is therefore convenient to 
study the ratio of the two contributions, 

R=°¥- ( 152 ) 

U ZM 

Up to small corrections this ratio is also independent of the flavour of the propagating 
states. The Higgs KK contribution can then be obtained by multiplying the zero-mode 
result by the corresponding R. It should be noted that not all topologies shown below 
actually contribute to in a specific RS model, either because the combination of SU(2) 
and U(l) group factors vanishes or because the model does not have Feynman rules that 
allow for the particular diagram to exist. Note that we do not separate contributions 
from wrong- and correct-chirality Higgs couplings. The numbers refer to the sum of both 
type of contributions, and are given by: 
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R ~ 0.27(0.01) 


R « 0.27(0.01) 




i ' i 


A 



R « 0.08(0.01) 



A « 0.08(0.01) 


7 


A « 0.15(0.05) 




A « 0.77(0.08) 


(153) 


Our error estimates are shown in parenthesis. It arises from the numerical integration 
error and an estimate for the extrapolation error to j3 = oo, since a numerically stable 
evaluation is possible only up to (3 ~ 200 — 300. The numerically most challenging 
diagrams are the ones where the photon is emitted from the Higgs, since they contain 
products of KK Higgs propagators. We also note that the KK Higgs contribution con¬ 
verges relatively slowly as [3 —> oo, if the diagrams involve an external Higgs attachment 
to a fermion line in the loop, as illustrated in Figure [2l] The typical scaling with powers 
of /3 in the different momentum regions discussed in the main text does not set in until 
/3 ~ 40. This behaviour agrees with observations made in 20 . 
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